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I.  INTRODUCTION 

The  weight  optimization  computer  program  OPTIM  II,  Reference  1, 
was  altered  to  make  its  input  format  compatible  with  the  widely  used 
NASTRAN  program.  The  resulting  new  program  is  called  OPTIM  III. 
Accomplishment  of  this  objective  required  the  modification  and  addition 
of  many  subroutines. 

Included  in  the  new  program,  OPTIM  III,  is  the  capability  of  analyzing 
and  optimizing  structures  which  contain  composite  membrane  plates.  The 
method  of  approach  adopted  for  analysis  of  composite  membrane  plates 
involves  the  stacking  of  triangle  and  quadrilateral  orthotropic  membrane 
plate  elements  and  is  discussed  in  Section  II.  Also  included  in  Section  II 
is  a  summary  of  the  three  types  of  failure  criteria  for  composite  materials 
in  a  state  of  plane  stress  which  are  available  in  OPTIM  III  for  use  in 
connection  with  the  optimization  process. 

The  theoretical  bases  for  the  triangle  and  quadrilateral  orthotropic 
membrane  plates  inserted  into  OPTIM  III  is  given  in  Section  III  which  also 
includes  a  summary  of  all  the  eight  finite  elements  available. 

All  NASTRAN  compatible  input  data  cards  are  described  in  Section  IV 
and  their  utilization  is  illustrated  by  a  sample  problem  in  Section  / . 

A  programmer's  manual  which  describes  the  computer  program 
OPTIM  III  in  detail  is  presented  in  the  appendix.  The  information 
presented  in  the  appendix  is  geared  specifically  to  the  programmer  and 
includes  the  general  program  logic,  external  file  structure  and  detailed 
descriptions  of  the  subroutines. 


OPTIMIZATION  OF  STRUCTURES  CONTAINING  COMPOSITE  PLATES 


A.  Genera  1 

The  membrane  triangle  and  quadrilateral  plate 
element  originally  included  in  OPTIMII  were  extended  to  handle 
composite  plate  constructions.  The  type  of  composite  struc¬ 
tures  accepted  by  the  optimization  program  is  limited  to  flat 
composite  plates  composed  of  crthotropic  layers  of  material. 
Each  layer  of  the  plate  composite  is  idealized  as  a  separate 
finite  membrane  plate  element  which  is  discussed  later  in 
this  section. 

When  modeling  a  composite  plate  consisting  of  many 
layers,  the  gridpoints  of  all  the  membrane  elements  or  layers 
of  the  plate  are  selected  to  be  the  same.  This  restriction 
ensures  that  each  of  the  layers  of  the  composite  plate  is 
subjected  to  the  same  deformation  state  and  is  in  conformance 
with  small  displacement  plate  theory  in  which  it  is  assumed 
that  straight  normals  to  the  plate  median  surface  before 
deformation  remain  straight  after  deformation.  To  ensu *e  that 
this  condition  is  satisfied,  the  layers  of  the  plate  composite 
must  be  oriented  symmetrically  with  respect  to  the  median 
surface  of  the  composite  plate. 

Various  failure  criteria  for  composite  materials 
were  investigated  and  of  these  three  types  of  criteria  were 
included  in  OPTIMIII.  These  three  types  of  criteria  which 
are  required  in  the  optimization  process  for  the  layers  in 
the  composite  plate  are  discussed  in  this  section.  (Reference 
6,  7,  8). 


B .  Analysis  with  Orthotropic  Membrane  Elements 

Each  layer  in  a  composite  membrane  plate  is  assumed 
to  be  composed  of  fibers  imbedded  in  a  matrix  material.  The 
fibers  are  assumed  to  be  oriented  in  such  a  way  as  to  result 
in  their  characterization  as  an  orthotropic  material.  A 
typical  orthotropic  layer  of  this  type  is  illustrated  in 
Figure  1.  Coordinates  Xm,  Ym  shown  in  the  figure  are  the 
material  axis  of  orthotropy  in  which  the  Xm  axis  is  oriented 
parallel  to  the  fibers. 

Two  types  of  orthotropic  membrane  plate  finite 
elements  are  available  in  the  optimization  program  and  they 
are , 

1.  Triangle  membrane  element 

2.  Quadrilateral  membrane  element 

These  two  finite  elements  are  displayed  in  Figure  2.  The  angle 
0  gives  the  orientation  of  the  material  axis  with  respect  to 
the  side  connecting  gridpoints  1  to  2.  Note  that  for  the 
triangle  the  material  axis  orientation  angle  0  is  measured 
from  side  1-2  to  the  material  axis  Xm  whereas  for  the  quadri¬ 
lateral  it  is  measured  from  the  material  axis  Xm  to  the  side 
1-2. 

The  relationship  between  the  elastic  stresses  and 
strains  for  an  orthotropic  material  in  plane  stress  referenced 
to  the  material  axis  of  orthotropy  shown  in  Figure  1  is  given  by 
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These  stress-strain  relations  were  used  in  the  development  of 
the  membrane  plate  elements  and  reduced  down  to  isotropic 
material  properties.  This  situation  makes  it  possible  to 
construct  a  composite  membrane  plate  composed  of  fiber- 
reinforced  orthotropic  materials  and  isotropic  metal i i c-type 
materials . 
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Figure  1.  Typical  fiber-reinforced  layer  in 
a  composite  membrane  plate 

The  stresses  are  computed  for  both  element  types 
with  respect  to  the  material  axis  of  orthotropy  and  are  used 
in  the  failure  criteria  formulation  discussed  in  the  next  section. 


Triangle 


b.  Quadrilateral 


Figure  2.  Orientation  of  membrane  elements  in 
fiber-reinforced  composite  layer 


C.  Criteria  of  Failure 

Three  types  of  criteria  of  failure  have  been  inserted  into  the 
optimization  program.  They  are 

a.  Modified  Von  Mises  or  energy  of  distortion  criterion 
of  failure 

b.  Maximum  Stress  Theory 

c.  Maximum  Strain  Theory 

In  ail  of  these  theories,  the  applied  stresses  c/y  an(*<yxy 
are  referred  to  the  material  axes  of  orthotropy. 

The  modified  distortion  energy  criterion  of  failure  is  given  by 


where  Fx  and  Fy  are  tensile  (compression)  failure  stresses  and  F^y  is  the 
shear  failure  stress.  These  failure  stresses  are  determined  for  simple 
load  test  conditions  with  respect  to  the  axis  of  orthotropy. 

The  maximum  stress  theory  postulates  that  failure  will  occur 
when  any  one  of  the  three  applied  stresses  is  equal  to  its  corresponding 
failure  stress.  Thus  failure  occurs  when  any  one  of  the  following  conditions 
is  satisfied 

<4  = 
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<*xy  ^xy 
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In  the  above  two  theories  tension  or  compression  failure  stresses  are 


selected  in  consistancy  with  the  sign  of  the  applied  stresses.  Failure 
stresses  Fx,  Fy  and  F^  represent  longitudinal,  transverse,  and  shear 
modes  of  failure. 

For  the  maximum  strain  theory  it  is  assumed  that  failure  is 
precipitated  when  any  one  of  the  applied  strain  components  referred  to  the 
material  axis  attains  its  limiting  failure  strain.  Thus  for  an  orthotropic 
material  the  conditions  of  failure  are 


where  £  ,  £  and  are  the  failure  strains.  Note  that  a  negative  applied 
x  y  xy 

strain  implies  a  compression  failure  strain.  The  condition  of  failure  in 
terms  of  stresses  are 
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Implementation  of  the  above  three  failure  criteria  requires 


a  knowledge  of  the  following  ten  material  failure  properties:  ±F  ,  ±F  , 

x  y 

Exy*  *^x»  ^Sy*  % xy' 
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Ill .  NEW  ORTHOTROPIC  PLATE  ELEMENTS 
A.  Gene  ra 1 

As  noted  in  Reference  1,  the  usefulness  of  an  optimization 
program  stems  from  the  accuracy  with  which  the  real  structure  can  be 
represented  by  the  finite  element  model.  In  the  OPTIM  1 1  program  there  are 
eight  elements  available  for  the  analysis  of  a  structure  fabricated  from 
isotropic  materials.  A  major  task  of  the  present  work  was  not  only  to 
put  OPTIM  II  into  a  form  compatible  with  NASTRAK  but  to  extend  its 
capability  to  include  composite  plate  structures.  In  this  latest 
version  of  OPTIM  n  there  are  eight  elements  as  shown  in  Figure  3  with 
their  program  names  and  they  are: 

1 .  Axial  Force  Member 

2.  Pure  Shear  and  Shear  Web  Elements 

3.  Triangular  Plate  in  Plane  Stress 

4.  Quadrilateral  Plate  in  Plane  Stress 

5.  Tubular  Beam  Element 

6.  Midpoint  Axial  Force  Member 

7.  Midpoint  Triangular  Plate  in  Plane  Stress 

8.  Midpoint  Quadrilateral  Plate  in  Plane  Stress 

The  first  three  elements  in  the  above  list  were  in  the  pilot 
optimization  program,  OPTIM  (Reference  2)  and  the  remaining  five  added 
when  OPTIM  II  was  written.  Explicit  formulations  for  the  element  stiffness 
and  other  matrices  are  given  in  Reference  2  for  the  first  five  elements 
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in  the  above  list.  Formulations  for  the  last  three  elements  are  given  in 


Reference  1.  The  triangular  and  quadrilateral  membrane  plate  elements, 
number  3  and  4,  were  extended  to  include  orthotropic  plate  properties 
and  the  details  of  this  extension  are  given  in  the  next  two  sections. 

B.  Triangular  Membrane  Plate 

The  triangular  plate  in  plane  stress  originally  contained  in 
the  OPTIM  II  program  was  extended  to  include  orthotropic  material 
properties.  The  location  of  the  material  axis  of  orthotropy  (X^  Ym) 
is  shown  in  Figure  3C<and  is  referred  to  the  local  reference  axes  (X,  Y) 
of  the  triangle  by  the  angle  0  .  A  description  of  the  procedure  used  to 
compute  the  stiffness  and  stress  matrices  is  given  below.  Theoretical  back¬ 
ground  for  the  procedure  is  given  in  References  2,  3  and  5. 

As  shown  in  Reference  2,  the  stiffness  matrix  referred  to  the 
local  axis  system  is  given  by 

[K]  =  t [B-1]1,  [A*]  [B"*]  (3.1) 

where 
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The  matrix  [A*]  is  given  by 

[A*]  =  [T]  [Ge]  [T]T 
where  the  transformation  matrix  [T]  is  .* 


(3.4) 
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For  isotropic 
given  by 
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l-M* 

[Ge]  =  ME 


rix  [Ge]  in  Equation  (3.4)  is 


(3.6) 


and  for  orthotropic  material; 

[O.J  =  [U]T  [Om]  [O] 


(3.7) 


where  the  rotation  transformation  matrix  [U]  of  strain  components  is  given  by 
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The  elements  cf  the  material  matrix  [Gm]  are  defined  in  terms  of  ortho¬ 
tropic  material  properties  as  follows: 
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The  membrane  stresses,  assumed  constant  in  the  element, 
are  given  with  reference  to  the  material  axis  of  orthotropy  by 


where 

[S]  =  [D]  [B'1] 

[D]  =  [TJ  [T]  [Ge]  [T]T  [B-1] 


(3.10) 

(3.  11) 

(3.  12) 

(3.13) 
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l.  Quadrilateral  Membrane  Plate 


The  quadrilateral  membrane  plate,  Figure  3d,  is  a  linear 
stress  variation  element  originally  developed  by  Turner^  for  isotropic 
materials  (9=0)  and  used  subsequently  in  many  analysis  and  optimization 
programs.  For  this  element,  the  design  variable  is  the  plate  thickness 
and  the  stress  variations  are  given  by 

tfx  =  al+a2^ 

( "5  -  a^+a^x 

Xxy  *  *5 

This  element  was  extended  to  include  orthotropic 
material  properties,  and  theoretical  details  of  the  extension  are 
given  in  Reference  3.  However,  the  element  matrices  developed  in 
Reference  3  were  referred  to  the  material  axis  and  to  make  them  compatible 
with  the  OPTIM  II  program^material  to  local  axis  transformation  had  to  be 


introduced. 


The  stiffness  matrix  for  the  quadrilateral  membrane  plate 


referred  to  the  local  element  axis  is  given  by 

[K]  =  t[T2]T  [B- 1  ]T  [A*]  [B_1]  [Tg] 


(3.14) 


The  matrix  [T2]  transforms  the  element  displacements  from  the  material 
axis  of  orthotropy  to  the  local  element  axis  and  is  given  by 
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(3.  17) 


Membrane  stresses  in  the  quadrilateral  membrane  plate  element 
are  given  by  Equation  3.  10  with  the  stress  matrix  defined  as 

[S]  =  [D]  [B-1]  [X2]  (3-18) 

where  matrix  [B]  is  given  by  Equation  3- 15 ,  transformation  matrix  [T2] 
by  Equation  3-17  and  the  [D]  matrix  by  the  following: 


[D]  = 


10  0  0 
0  0  10  0 


Ym  0  0 


P  0  0  xm  0 


0  0 
0  0 


(3.  19) 


where  Xm,  Ym  are  referred  to  the  material  axis. 
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Stresses  used  in  the  failure  criteria  for  the  quadrilateral 


IV.  PROGRAM  INPUT  DESCRIPTION 


A.  Input  Overview 

This  section  describes  the  input  data  required  to  execute 
the  NAS TRAN  compatible  optimization  program  OPTIM  III. 

OPTIMIIIwas  developed  such  that  most  of  the  input  data 
cards  are  in  NASTRAN  format.  This  input  consists  of  a  deck 
beginning  with  "BEGIN  BULK"  and  ending  with  "ENDQATA".  All 
data  cards  are  optional  except  those  defining  controls,  loads, 
grid,  and  boundary  conditions.  Elements  must  also  be  defined. 

The  NASTRAN  input  feature  of  submitting  data  in  any 
order,  either  left  or  right  adjusted,  is  preserved.  The 
GRID  cards,  SPC  cards,  FORCE,  MOMENT,  and  MATT,  MAT2  cards 
are  similar  to  NASTRAN  input.  Figure  3  shows  the  OPTIM  III 
element  library  and  the  appropriate  CrTIM  input  code  for  each 
element.  The  following  is  a  summary  of  OPTIM  III  Input  Cards: 


B.  Definitions  of  Input  Cards 


BEGIN  Bulk 
♦Buck 
♦Buck! 
CONROD 
CQUAM8 
CQDMEM1 
CRDMID 
CROD 
CSNEAR 
CTRIM6 

CTRMEM 

CTUBEAM 

CUES 

FORCE 

♦GCON 

♦GPCEQN 

GRID 

♦ICON 

♦LINKS 

MAT! 

MAT2 

MOMENT 

♦OPDVIR 

OPLOADS 

OPTIM 

PQUAM8 

PQDMEM1 

PRDMID 


First  Card  of  Input 

Define  elements  which  require  buckling 
Define  set  of  elements  which  require  bucklino 
AXIAL  FORCE  MEMBER  Property  and  Connection 
8  Noded  Midpoint  Quadrilateral  Connection 
Quadrilateral  Membrane  Connection  Card 
Midpoint  AXIAL  Force  Member  Connection 
Axial  Force  Element  Connection 
Shear  Panel  Element  Connection 
Triangular  6  node  midpoint  membrane  element 
connection 

Triangular  membrane  element  connection 
Tubular  beam  element  connection 
Shear  Web  element  connection 
Static  load  at  grid  point 
Generalized  Constraint  Components 
Grid  Point  Constraint  Equations 
Grid  Point  Coordinates 
Individual  Constraints 
Linked  group  element  members 
Isotropic  material  properties 
Anisotropic  material  properties 
Static  moment  at  a  gridpoint 
Selected  element  design  variables 
Loads  for  optimization 
Optim  control  parameters 

Property  card  for  8  node  midpoint  quadrilateral 
Property  card  for  quadrilateral  membrane 
Property  card  for  midpoint  axial  force  membrane 


♦These  cards  are  optional  and  need  to  be  input  only  if  necessary 
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PRO"  Property  card  for  axial  force  member 

PSml.hR  Property  card  for  shear  panel 

PTRIM6  Property  card  for  triangular  6  node  midpoint 

PTRMEM  Property  card  for  triangular  membrane 

PTUBEAM  Property  card  for  tubular  beam 

PWEB  Property  card  for  shear  web 

SPC  Single  point  constraint 

SPC1  Sets  of  single  point  constraints 

SPCI  THRU  Sets  of  single  point  constraints 

TITLE  Title  card  information 

ENDDATA  End  of  Data  Deck 


C.  Input  Format 

In  this  section  detail  OPT I M  III  input  card  descriptions 
are  presented  in  a  manner  similar  to  the  NASTRAN  User's  Manual. 
That  is,  each  data  card  is  described  individually  in  alphabet¬ 
ical  order. 


Input  Data  Card  BEGIN  Bulk 

Description :  First  Card  of  Input 

I  23456789  10 

BEGIN  Bulk 


This  card  is 
in  a  NASTRAN 


identical  to  the  format  of  the 
bulk  data  deck. 


first  card  used 


Input  Data  Card  BUCK 


Buckling  Elements 


Descri pti on :  Defines  elements  which  require  buckling 

*  ”  calculations 


Format  and  Example: 


1 

8UCK 

2 

SID 

3 

EL] 

4 

el2 

5 

el3 

6 

el4 

7 

el5 

8 

el6 

9 

ELy 

BUCK 

1 

2 

3 

7 

8 

9 

10 

12 

Field 

SID 

EL-| ,  etc. 


Identification  number  of  the  buckling  set 
Element  number 


Input  Data  Card  BUCK! 


Buckling  Elements 


Description:  Definessets  of  element  which  require  buckling 

cal cul ati ons 


Format 

and  Example: 

1 

2 

3 

BUCK! 

SID 

EL 

BUCK! 

1 

1 

Field 

SID 

Identification  i 

ELl 

Element  number 

ELn 

Element  number 

"  THRU" 
THRU 


the  buckling  set 


l  f 
-  !\ 

!  ■ 

i 


'i 


BULK  DATA  DECK 


Input  Data  Card  C0NR0D  AXIAL  FORCE  MEMBER  Property  and 

Connection 


Descri pti on :  Defines  a  rod  element  of  the  structural  model 

without  reference  to  a  property  card. 


Format  and  Example: 


12  3  4 

C0NR0D  EID  G1  G2 

C0NR0D  2  16  17 


5  6  7 

MID  A 

23  2.69 


8  9 


10 


Field  Contents 

EID  Unique  element  identification  number  (lnteger>0) 

G1 ,  G2  Grid  point  identification  numbers  of  connection 
points  (lnteger>0;  G1  f  G2) 

MID  Material  identification  number  (lnteger>0) 

A  Area  of  rod  (Real) 


Remarks :  1.  Element  identification  numbers  must  be  unique 

with  respect  to  al 1  other  element  identification 
numbe  rs . 

2.  For  structural  problems,  C0NR0D  cards  may  only 
reference  MAT!  material  cards. 
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Input  Data  Card  CQUAH8  Midpoint  Quadrilateral  Element  Connection 


Description;  Defines  a  midpoint  quadrilateral  membrane  element 
(QUAM8)  of  the  structure. 


Format  and  Example: 


1 

2 

3  4 

5 

6 

7  8 

9 

1 

CQUAM8 

EID 

PID  G 1 

G2 

G3 

G4  G5 

G6 

+a 

CQUAM8 

10 

3  1 

10 

11 

8  9 

20 

S 

1 

2 

3  4 

5 

6 

+abc 

G7 

G8 

+  ST 

13 

16 

Field 

Contents 

EID 

Element  identification 

number 

(Integer  >0) 

PID 

Identi fi cation 

number 

of  a  PQUAM8  property 

card 

G1  through  Gridpoint  identification  numbers  of  connection  points 
G8  (Integer)O) 


Remarks :  1.  Element  identification  numbers  must  be  unique 

with  respect  to  al 1  other  element  identification 
numbers. 

2.  Gridpoints  Cl  through  G8  must  be  ordered  as  shown 
in  the  sketcr  of  the  element. 


Input  Data  Card  CQDMEM1  Quadrilateral  Element  Connection 


Descri pti on :  Defines  a  quadrilateral  membrane  element  (QDMEM1) 
of  the  structural  model. 


Format  and  Example: 

1  2  3 

CQDMEM1  El D  PID 

CQDMEM1  72  13 


4  5  6 

G1  G2  G3 

13  14  15 


7  8  9 

G4  TH 

16  29.2 


Field 


Contents 


EID  Element  identification  number  (Integer^O) 

PID  Identification  number  of  a  PQDMEM1  property  card 

(Integer  >0) 


G1,G2,  Grid  point  identification  numbers  of  connection  points 
G3,G4  (Integer  >0) ;  G1  f  62  f  G3  f  G4) 


TH  Material  property  orientation  angle  in  degrees  (Real). 

The  sketch  below  gives  the  sign  convention  for  TH. 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  al  1  other  element  identification  numbers 

2.  Gridpoints  G1  through  G4  must  be  ordered  con¬ 
secutively  around  the  perimeter  of  the  element. 


3.  All  interior  angles  must  be  less  than  180  degrees. 


Input  Data  Card  C  RDM  I D  Midpoint  Axial  Force  Element  Connection 


Description :  Defines  a  midpoint  axial  force  element  RDMID  of 
the  structural  model. 


Format  and  Example: 

1  2  3 

CRDMID  EID  PID 

C RDM  ID  72  13 


4  5 

G1  G2 

13  14 


6  7  8  9 

G3 

15 


Field 


Contents 


EID 


Element  identification  number  (lnteger>0) 


PID  Identification  number  of  a  PRDMID  property  card 

(Integer  ^  0) 

G1,G2,G3  Gridpoint  identification  numbers  of  connection 

points  (Integer^  0);  G1  f  G2  f  G3 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  al  1  other  element  identification  numbers 

2.  Gridpoints  G1  through  G3  must  be  ordered  as  shown. 


Input  Data  Card  CR0D 


Axial  Force  Element  Connection 


Description:  Defines  a  tension-compression  element  (R0D)  of 

the  structural  model 


Format  and  Example: 


1 

2  3 

4 

5 

6  7 

8 

9  10 

CR0D 

EID  PID 

Gl 

G2 

EID  PID 

Gl 

G2 

CR0D 

12  13 

21 

23 

3  12 

24 

5 

Field 

Contents 

EID 

Element 

i denti f i cati on 

number  (Integer  >  0) 

PID 

Identifi cation 

number 

of  a  PR0D 

property 

card 

Gl,  G2  Grid  point  identification  numbers  of  connection  points 
(Integer  >  0 ;  Gl  f  G2 ) 


1.  Element  identification  numbers  must  be  unique  with 
respect  to  al 1  other  elem.  t  identification  numbers 

2.  See  C0NR0D  for  alternative  method  of  rod  definition 

3.  One  or  two  R0D  elements  may  be  defined  on  a  single 
card. 


Remarks 


•1 


Input  Data  Card  CSHEAR  Shear  Panel  Element  Connection 

Descri pti on :  Defines  a  shear  panel  element  (SHEAR)  of  the 
structural  model. 


Format  and  Examt 


1  2  3 

CSHEAR  EID  PID 

CSHEAR  3  6 


4 

G1 

1 


5  6  7  8 

G2  G3  G4 

5  3  7 


9  10 


Field 

EID 

PID 

Gl,  G2 , 
G3,  G4 


Contents 

Element  identification  number  (Integer^  0) 

Identification  number  of  a  PSHEAR  property  card 

Grid  point  identification  numbers  of  connection  points 
( Integer>  0;  Gl  f  G2  f  G3  f  G4) 


Remarks :  1.  Element  identification  numbers  must  be  unique 

with  respect  to  all  other  element  identification 
numbers . 

2.  Grid  points  Gl  thru  G4  must  be  ordered  consecutively 
around  the  perimeter  of  the  element. 


-  f 
1 


Input  Data  Card  CTRIM6  Midpoint  Triangular  Membrane  Plate 

Element  Connection 


Description:  Defines  a  triangular  membrane  element  ( TRIM6 )  of 
the  structural  model. 


Format  and  Example: 


1  2  3  4  5 

CTRIM6  EID  PID  G1  G2 

CTRIM6  220  666  100  110 

+abc 

+BS  90.0 


6  7  8 

G3  G4  G5 

120  210  220 


9 

G6  h 
320 


Field 


Contents 


EID 

PID 

G1-,G2,G3, 

G4,G5,G6 


Element  identification  number  (lnteger>0) 

Identification  number  of  PTRIM6  property  card. 

Gridpoint  i denti fi cati on . numbers  of  connection 
points  (Integers >  0) ;  G1  G2  *  G3  G4  ?  G5  G6). 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  a  1 1  other  element  identification  numbers 

2.  Interior  angles  must  be  less  than  180°. 

3.  The  gridpoints  must  be  ordered  consecutively  around 
the  perimeter  in  a  counter  clockwise  direction  and 
starting  at  a  vertex. 


4.  The  continuation  card  must  be  present. 


5.  Gridpoints  G2,  G4  and  G6  are  assumed  to  lie  at 
the  midpoints  of  the  sides. 


Input  Data  Card  CTRMEM 


Triangular  Element  Connection 


Descri pti on : 

Defines  a  triangular 

membrane 

element  (TRMEM)  of 

the  structural 

model . 

Format  and  Example: 

1  2 

3  4 

5 

6 

7  8  9 

CTRMEM  EID 

PID  G1 

G2 

G3 

TH 

CTRMEM  16 

2  12 

1 

3 

16.3 

■i 


i  j 


i 1 


Field 

EID 

PID 

G1 , G2  ,G3 
TH 


Contents 


Element  identification  number  (Integer  ^0) 

Identification  number  of  a  PTRMEM  property  card 

Gridpoint  identification  numbers  of  connection  points 
(Integer  >0;  G1  t  G2  f  G3) 

Material  property  orientation  angle  in  degrees 
(Real)  -  The  sketch  below  gives  the  sign  convention 
for  TH. 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  all  other  element  identification  numbers. 


,41 

4  ; 

‘4  s 
J  ' 

FT  i 

i?  { 


i 


Input  Data  Card  CTUBEAM 


Tubular  Beam  Element  Connection 


Defines  a  tension-compression-torsion  element 
(TUBE)  of  the  structural  model. 


Format  and  Example: 


1  2  3  4  5  6 

CTUBEAM  El D  PID  G1  G2  G3 

CTUBEAM  12  13  21  23  3 


Field  Contents 

EID  Element  identification  number  (lnteger>0) 

PID  Identification  number  of  a  PTUBEAM  property  card 

(Default  is  EID)  (lnteger,>0) 

G1 ,  G2,  G3  Gridpoint  identification  numbers  of  connection  points 
(lnteger>0;  G1  y*  G2  ^  G3) 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  al 1  other  element  identification  numbers 

2.  See  Figure  3e. 


Input  Data  Card  CWEB  Shear  Web  Element  Connection 

Descri ption :  Defines  a  2  node  symmetric  Web  element  (Web) 

of  the  structural  model 

Format  and  Example: 

1  2  3  4  5  6  7  8  9 

CWEB  EID  PID  G1  G2 

CWEB  4213 

Field 

EID  Element  identification  number 

PID  Identification  of  PWEB  property  card 

G1 ,  G2  Gridpoint  identification  numbers 


Remarks:  1.  See  Figure  3b, 


Input  Data  Card  GCON  Generalized  Constraint  Components 


Descri pti on :  Defines  the  translation  and  rotation  components 

of  each  gridpoint  which  are  active  degrees  of 
freedom  in  a  generalized  constraint  equation 
defined  by  GPCEQN  cards 


Format  and  Exami 


1234  5  6789  ID 

gcon  sid  gp  luiviw  ±exi9yte2 
GCON  14  1  +XA-X  A+X  A 


Field 

SID  Equation  identification 

GP  Gridpoint  number 

+,-  Algebraic  sign  of  component  coefficient 

x  The  x  mark  indicates  that  the  comoonent  exists. 

A  A  triangle  indicates  that  the  component  does  not 

exist. 

Remarks :  The  SID  numbers  on  these  cards  appear  in  the  GPCEQN  cards 

which  set  up  the  constraint  equation.  See  remarks 
under  GPCEQN  for  further  explanation. 


Input  Data  Card  GPCEQN  Generalized  Constraint  Equation 


Description:  Defines  the  upper  and  lower  limits  of  a  generalized 

constraint  equation.  The  existance  of  each  GPCEQN 
card  is  necessary  for  each  constraint  equation. 


Format  and  Example: 

1234567891 
GPCEQN  SID  GCS  Upper  Lower 

LimityL  LimityL 

GPCEQN  1  14  4.0  1.0 


Field 

SID  Equation  identification 

GCS  Set  identification  of  GCON  cards  which  define  the 

components  in  the  equation. 

UL  Upper  limit  of  constraint  equation 

LL  Lower  limit  of  constraint  equation 


Remarks :  1.  Suppose  the  user  wants  to  define  the  following 

generalized  constraint  equations: 


1 .0<  -U]  +  V2  +  0xl  £4.0 

1.0  <  W2  -0x4<6.0 

The  input  cards  would  be  as  follows: 


1 

2 

3 

4 

_ 5 _ 

GPCEQN 

1 

14 

4.0  _  ... 

1.0- 

GPCEQN 

2 

12 

fi.O 

1 .0  . 

GCOfi 

14 

1 

-X 

+x 

GCON 

14 

2 

+x 

GCCN 

12 

2 

+x 

GCON 

12 

4 

-X 

Equation  #1 
Equation  #2 


-«1.  «xi 

V2 

Wo 

x4 


2.  A  generalized  constraint  consists  of  the  algebraic 
combination  of  the  selected  degrees  of  freedom  of 
all  participating  gridpoints  and  is  maintained 
within  the  upper  and  lower  limits  defined  for 
the  particular  group. 


Input  Data  Card  GRID 


Grid  Point 


Descri pti on :  Defines  the  location  of  a  geometric  gridpoint  of 

the  structural  model  and  its  permanent  single¬ 
point  constraints. 


Format  and  Example: 


1  2  3 

GRID  ID 

GRID  2 


4  5  6 

XI  X2  X3 

1.0  2.0  3.0 


7 


8 

PS 

316 


9 


10 


Contents 

Gridpoint  identification  number  (0<Integer<999999 ) 
Location  of  the  grid  point. 

Permanent  single-point  constraints  associated  with  grid- 
point  {any  of  the  digits  1-6  with  no  imbedded  blanks) 
(Integer>  0  or  blank) 


Remarks :  1.  The  coordinate  system  defined  on  all  GRID  cards  is 

called  the  Global  Coordinate  System.  All  degrees- 
of-freedom,  constraints,  and  solution  vectors  are 
expressed  in  the  Global  Coordinate  System. 


Field 

ID 

XI ,X2,X3 
PS 


i 

( 

r 

i 

5 

I 


Input  Data  Card  I  CON  Individual  Constraints 


Des cri pti on :  Defines  the  constraint  limits  (both  lower  and 

upper)  and  gridpoint  comoonents  which  require 
individual  constraints. 


Format  and  Example: 


}  1  2 

3  4 

5 

6 

7 

8 

9 

)  ICON  Components 

Lower  Upper 

Si 

g2 

S3 

S4 

g5 

l  ICON  146 

• 
r 

l 

i 

1.0  3.0 

1 

2 

3 

6 

\  Field 

1 

f  Components 

1-6  signifies 

component 

numbers 

Lower  Lower  limit 


10 


Upper  Upper  limit 

G i  > G £ » G 3 »  etc:*  Gridpoint  numbers.  Upto  5  gridpoint  numbers 

may  appear  on  one  ICON  card. 


'4 


1 

*r 


f 


Remarks : 


Each  ICON  card  contains  up  to  5  gridpoints.  When  more  than 
5  are  desired,  list  them  on  separate  ICON  cards. 

For  example,  let  there  be  the  following  individual  constraints 
U i ,  Wi »  ^2  Sweater  than  1.0  and  less  than  4.0. 

V],  Vg,  V3  greater  than  1.0  and  less  than  8.0. 


The  required  ICON 

cards 

would  be 

* 

• 

1  2 

3 

4 

5 

6 

7 

ICON  13 

1.0 

4.0 

1 

2 

I  CON  2 

1.0 

8.0 

1 

2 

3 

3.  Rotational 
element  is 

1  i  mi ts 
used. 

are  indicated 

only 

when 

3ft 


^ B 


Input  Data  Card  LINKS 


Descri pti on :  Defines  element  numbers  which  are  in  a  "linked" 

group  either  for  purposes  of  specifying  area  ratios 
or  for  buckling  definitions. 


Format  and 

Exampl e : 

1  2 
LINKS 

LINKS 

3  4  5  6 

BELEM  NG  BUCK  N-. 

5  3  YES  1 

7  8 

N2  N3 

2  3 

9 

10 

Field 

Contents 

BELEM 

Base  element  number  entered  in 

integer  format 

NG 

Number  of  elements  in  the  group.  The  total  number 
elements  including  the  base  element  is  entered. 

of 

BUCK 

Buckling  control.  When  "yes" 
the  grouping  defined  is  to  be 
and  n£  linking  between  design 

is  in  this 
considered 
variables 

field,  then 
for  buckling 
exist. 

N  t  ,  N  2  » N  3 

Node  points  defining  rectangular  area. 

Remarks : 

1.  This  data  is  optimal. 

2.  The  base  element  is  the  lowest  numbered  element 

in  the  group  and  it  must  also  be  the  element  with 
the  smallest  minimum  allowable  element  design 
variable.  All  elements  within  a  group  are 
directly  related  to  the  base  element  through 
the  area  ratios  specified  on  the  element  property 
cards.  In  other  words  the  minimum  allowable 
area  is  entered  for  the  base  element.  For 
other  members  in  the  group  the  ratio  (>1)  of 
member  area  to  base  element  area  is  entered 
instead  of  an  allowable.  This  ratio  is  main¬ 
tained  during  the  solution. 

3.  Tne  example  above  indicates  that  elements  5,  6 
and  7  belong  to  a  linked  group  in  which  element  5 
is  the  base  element. 

4.  The  buckling  option  has  been  included  in  order 
that  a  User  might  specify  a  simply  supported 
area  that  is  comprised  of  more  than  one  element. 


( conti nued) 


The  three  node  point  numbers  are  entered  as  fixed 
point  numbers  in  order  to  define  the  rectangular 
simply  supported  region.  The  sketch  illustrates 
the  three  points  and  their  relationship  to  the 
contained  groups  of  elements. 


Simply  supported"^  ^ 
perimeter 


Typical  quadrilateral 
element 


Input  Data  Card  M0MENT  Static  Moment 


Descri pti on :  Defines  a  static  moment  at  a  grid  point  by 

specifying  a  vector. 


Format  and  Example: 

1  2  2 

M0MENT  SID  G 

M0MENT  2  5 


M 

2.9 


6 

rn 

o.o 


7  8  9 

M2  N3 

1.0  0.0 


Field 

SID 

G 

M 

N1 , N  2 , N  3 


Contents 

Load  set  identification  number  (lnteger>0) 
Gridpoint  identification  number  (lnteger>0) 
Scale  factor  (Real) 
components  of  Vector 


Remarks : 


1.  The  static  moment  applied  to  gridpoint  G  is  given 
by 

m"  =  M •  (Ml , N 2 , N 3 ) 

2.  Load  set  is  selected  on  the  OPLOADS  card. 
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Input  Data  Card  HAH  Material  Property  Definition 

Descri pti on :  Defines  the  material  properties  for  linear, 

temperature-independent,  isotropic  materials. 


Format  and  Example: 


12  3  4 

MAT1  MID  E  G 

MAT!  17  3. +7  1.9+7 


5 

MU 


+abc  SL  SU 
+BC  20. +4  15. +4 


6  7 

RH0 
4.28 


8 


9  10 

+abc 
ABC 


Field 

MID 

E 

G 

NU 

RH0 

SL,  SU, 


Contents 

Material  identification  number  (Integer^O) 
Young's  modulus  (Real>0.0  or  blank) 

Shear  modulus  (Real>0.0  or  blank) 

Poisson's  ratio  (-1 . 0<  Real  <  0. 5  or  blank) 
Mass  density  (Real) 

Lower  stress  limit,  upper  stress  limit 


Remarks :  1.  One  of  E  or  G  must  be  positive  (i.e.,  either  E>0.0 

or  G>0.0  or  both  E  and  G  may  be>0.0). 

2.  If  any  one  of  E,  G  or  NU  is  blank,  it  will  be 
computed  to  satisfy  the  identify  E  -  2(1+NU)G; 
otherwise,  values  supplied  by  the  user  will  be 
used. 

3.  The  material  identification  number  must  be  unique 
for  all  MAT1  and  MAT2  cards. 
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Input  Data  Card  MAT2 


Material  Property  Definition 


Description:  Defines  the  material  properties  for  linear, 
anisotropic  materials. 


Format 

and  Example: 

(consists  of 

2  cards) 

1 

2  3 

4 

5 

6  7 

8 

9 

MAT2 

MID  Gil 

G 1 2 

G13 

G22  G23 

G33 

RHO 

MAT2 

13  6.2+3 

6.2+3 

5.1  +  3 

0.056 

+  abc 

SL 

SU 

+BC 

20. +5 

Field  Contents 

MID  Material  identification  number  ( Integer  ,>0) 

Gij  The  material  property  matrix  (Real) 

RH0  Mass  density  (Real) 

SL ,SU  Lower  stress  limit,  upper  stress  limit 


Remarks :  1.  The  material  identification  numbers  must  be  unique 

for  al 1  MAT! ,  MAT2  cards. 

2.  The  convention  for  the  G-jj  in  fields  3  through  8 
is  represented  by  the  matrix  relationship. 


3. 


Gn 

G1  2 

G1 3 

*') 

g12 

g22 

g23 

\  £2 

G13 

S23 

g33. 

lri2J 

Only  TRMEM  and  CQDMEM  elements  may  use  MAT2  cards. 
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Input  Data  Card  OPDVIR 


Descri pti on : 


Defines  design  variables  for  selected  elements, 
These  design  variables  are  used  as  starting 
guesses  in  an  optimization  run  or  else  as 
specified  values  in  a  status  run. 


Format  and  Example: 


OPDVIR 

OPDVIR 


3 

DVARi 

.1 


5  6 

DVARp  E3 
.6  4 


7  8  9 

DVAR3  E4  DVAR4 

.50  7  .23 


Field 


El ,E2. . .etc. 


DVAR-|  ,  DVAR2, 
etc . 


Element  numbers  for  which  design  variables  are 
i nput. 

Design  variable  starting  guesses  for  elements 
1 ,  2,  etc. 


Remarks :  1 . 


Up  to  4  elements  may  be  defined  on  one  OPDVIR 
card.  More  elements  may  be  defined  by  successive 
OPDVIR  cards. 


The  use  of  these  values  is  dependent  on  the 
calculation  control  on  the  OPTIM  input  card  - 
field  #9,  OPT. 

a.  If  OPT  =  GOPT  the  program  will  optimi ze 
using  the  design  variable  minimums  input 
in  this  section  as  a  starting  guess. 

b.  If  OPT  =  GSTA,  the  program  will  perform  a 
stati cs  solution  using  the  design  variables 
i nput  i n  this  section. 


I r  p ; t  Data  Card  OPLOADS 


Descri pti on :  Defines  the  character  of  the  loads  for  optimization 

i.e.,  symmetric,  nonsymmetric ,  antisymmetric. 


Format  and  Example: 


OPLOADS  LOAD  NSYM  SYA  F, 
OPLOADS  1  ANTI  5 


Mi  F2  M2 


Field 

LOAD 

NSYM 

SYA 


Fi 

Mi 


Contents 

LOAD  number  identification  (any  number) 

Number  of  nodes  on  symmetric  plane 

SYMM  -  Symmetric  loads 

ANTI  -  Antisymmetric  loads 

NON  -  Nonsymmetric  loads 

SID  number  of  FORCE  cards 

SID  number  of  MOMENT  cards 


F2  required  only  SID  number  of  FORCE  cards 

M2  for  nonsymmetric  SID  number  of  MOMENT  cards 


Remarks :  1.  If  Fi  is  present,  i  force  cards  are  required. 

J-r  M,;  is  present,  i  moment  cards  are  required. 


2. 


Each  external  load  requires  an  OPLOADS  entry. 


Input  Data  Card  OPTIH 

Description :  Defines  various  control  parameters  for  the  OPTIM 
program. 


Format  and  Example: 

1  2  3  4  5  6 

OPTIH  PRI  PRS  PRE  BU  MAX 

OPTIM  NO  YES  YES  NONE  10 


7  8  9 

CONVD  CONVW  OPT 

.01  100.0  OPT 


10 


Field  Contents 

PRI  NO/YES  Print  displacements  every  iteration 

PRS  NO/YES  Print  stresses  every  iteration 

PRE  NO/YES  Print  element  design  variables  every  iteration 

BU  NONE/ALL/SEL  Buckling  Anaysis  SEL=SID  of  BUCK 
MAX  Maximum  number  iterations 
CONVD  Convergence  criteria  for  design  number 
CONVW  Convergence  criteria  for  weight  perturbations 
OPT  OPT  Optimize  using  design  variable  on  property  card 
GOPT  Optimize  using  OPDVIR  starting  guess 
GSTA  Statics  run  using  OPDVIR 
STA  Statics  run  using  property  card 

Remarks :  1.  CONVD  -  Convergence  Criteria  (Design  Variable).  If 
two  successive  iterations  of  element  design  variables 
i.e.,  thickness  or  area,  meet  the  criteria  I  An'^n+1 

then  the  iteration  is  stopped.  C-j  is  entered  here.  ^n  ' 

2.  CONVW  -  Convergence  Criteria  (Weight  Perturbation). 

If  two  successive  perturbations  of  the  structural  weight  meet 
this  criteria,  then  the  iteration  is  stopped. 

Wj  c2 

“n 

Cj  is  entered  here. 
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Input  Data  Card  PQUAM8  MIDPOINT  QUADRILATERAL  MEMBRANE  PROPERTY 
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Input  Data  Card  PRDMID 


Descri pti on : 


Used  to  define  the  properties  of  a  midpoint 
axial  force  member.  Referenced  by  the  CRDMID  card. 


Format  and  Example: 


123456739  10 

PRDMID  PID  MID  A  PID  MID  A 

PRDMID  24  2  .49  25  3  .60 


Field 

PID  Property  identification  number  (integer>0) 
MID  Material  identification  number  (integer>0) 
A  Cross  sectional  area  (minimum) 


Remarks : 

1.  All  PRDMID  must  have  uniaue  property  identification  numbers. 

2.  One  or  two  midpoint  axial  force  member  properties  may  be 
defined  on  a  single  card. 


Inout  Data  Card  PR0D  Rod  Property 


Descri pti on :  Defines  the  properties  of  a  rod  which  is 

referenced  by  the  CR0D  card 


Format  and  Example 


! 

t 

> 

I 

* 


1  2  3 

PR0O  PID  MID 
PR0D  17  23 


4 

A 

42.6 


6  7  8  9  10 


Field  Contents 

PID  Property  identification  number  (Integer  >0) 

MID  Material  identification  number  (Integer^O) 

A  Area  of  rod  ( Real ) 


Remarks :  1  PR0D  cards  must  all  have  unique  property 

identification  numbers. 

2.  For  structural  problems,  PR0D  cards  may  only 
reference  MAT!  material  cards. 


J 
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. I 


Input  Data  Card  PSHEAR 


Shear  Panel  Property 


Descri pti on :  Defines  the  elastic  properties  of  a  shear  panel. 


Referenced  by 

the  CSHEAR 

card . 

Format 

and  Example: 

1 

2 

3 

4 

5  6 

7 

8 

9 

PSHEAR 

PID 

MID 

T 

PID 

MID 

T 

PSHEAR 

13 

2 

4.9 

14 

6 

4.9 

Field  Contents 

PID  Property  identification  number  (lnteger>0) 

MID  Material  identification  number  (Integer  >0) 

T  Minimum  thickness  of  shear  panel  (Real  f  0.0) 


Remarks :  1.  All  PSHEAR  cards  must  have  unique  identification 

numbers . 


2.  PSHEAR  cards  may  only  reference  MAT1  material  cards. 

3.  One  or  two  shear  panel  properties  may  be  defined 
on  a  single  card. 


Input  Data  Card 


PTRIM6 


Linear  Strain  Midpoint  Trianaular 
Membrane  Property 


Descri pti on :  Defines  the  properties  of  a  midpoint  triangular 

membrane  element.  Referenced  by  the  CTRIM6  card. 
No  bending  properties  are  included. 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

PTRIM6 

PID 

MID 

T1 

PTRIM6 

666 

999 

1.17 

Field 

PID 

Property 

i denti f i cati on 

number 

( Integer> 

0) 

MID 

Material 

i denti f i cation 

number 

(Integer^ 

0) 

T1 

Membrane 

thickness  of  the  element  (Real) 

t 


Remarks : 


1.  All  PTRIM6  cards  must  have  unique  property 
identification  numbers. 

2.  PTRIM6  cards  may  only  reference  MAT!. 


Input  Data  Card  PTRMEM 


Triangular  Membrane  Property 


Descri  pti  on :  Used  to  define  the  properties  of  a  triangular 

membrane  element.  Referenced  by  the  CTRMEM  card 
No  bending  properties  are  included. 


Format  and  Example: 


12  3  4 

PTRMEM  PID  MID  T 

PTRMEM  17  23  4.25 


Field  Contents 

PID  Property  identification  number  (lnteger,>0) 

MID  Material  identification  number  (lnteger,>0) 

T  Membrane  thickness  (Real>0;0) 


Remarks :  1.  All  PTRMEM  cards  must  have  unique  property 

identification  numbers. 

2.  One  or  two  triangular  membrane  properties  may  be 
defined  on  a  single  card. 


Input  Data  Card  PTUBEAM 


Tubular  Beam  Property 


Description:  Defines  the  properties  of  tubula“  beam  element. 

Referenced  by  the  CTUBE  card. 


Format  and  Example: 


1  2  3 

PTUBEAM  PID  MID 

PTUBEAM  2  6 


4  5  6 

TT  RAD 

6.29  0.25 


7  S  9 


Field 

PID 

MID 

TT 

RAD 


Contents 

Property  identification  number  (Integer  >0) 
Material  identification  number  (lnteger>0) 
Minimum  tube  thickness 
Radius  of  tube  element 


Remarks :  1.  PTUBE  cards  must  all  have  unique  property 

identification  numbers. 


2.  PTUBE  cards  may  only  reference  MAT!  material  cards 


Input  Data  Card  PWEB  WEB  Property 

Description:  Used  to  define  the  properties  of  2  node  WEB. 

Referenced  by  the  CWEB. 


Format  and  Example: 


1  23456789  10 

PWEB  PID  MID  T  PID  MID  T 

PWEB  235  2  0.5 


Contents 

Property  identification  number  (Integer^  0) 
Material  identification  number  (Integer  >0) 
Thickness  of  the  shear  web  (Real  >0.0) 

Remarks :  1.  All  PWEB  cards  must  have  unique  property 

identification  numbers. 

2.  One  or  two  WEB  properties  may  be  defined  r>n  a 
single  card. 


Field 

PID 

MID 

T 
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Input  Data  Card  SPC 


Single-Point  Constraint 


Descri pti on :  Defines  sets  of  single-point  constraints 


Format  and  Example: 


1  2  3 

SPC  SID  G 

SPC  2  32 


4  5 

C 

436 


6 

G 

5 


7  8  9  1C 

C 
1 


£ je 1 d  Contents 

SID  Identifi cation  number  of  single-point  constraint  set 

( Integer >  0) 

G  Grid  point  identification  number  ( Integer >0 ) 

C  Component  number  (Any  unique  combination  of  the  digits 

1-6  (with  no  imbedded  blanks)  when  point  identification 
numbers  are  gridpoints). 


Remarks : 

1 . 

Single-point  constraint  sets  must  be  present  in 
the  input. 

2. 

From  one  to  twelve  single-point  constraints  may 
be  defined  on  a  single  card. 

3. 

SPC  degrees  of  freedom  may  be  redundantly 
specified  as  permanent  constraints  on  the 

GRID  card. 

4.  The  SID  number  should  be  the  same  on  all  SPC 


Input  Oata  Card  S P C 1  Single-Point  Constraint 


Oescri pti on :  Defines  sets  of  sinale-point  constraints. 


Format  and  Example: 

123456789  10 

S  PC  1  SID  C  G1  G2  G3  G4  G5  G6 

SPC1  321  3  10  9  6  5 

Alternate  Form 

SPC1  SID  C  GID1  “THRU"  G I D2 
SPC1  313  12456  6  THRU  32 


Contents 

Identification  number  of  single-point  constraint  set 
( Integer  >0) 

Component  number  (Any  unique  combination  of  the  digits 
1-6  (with  no  imbedded  blanks)  when  point  identification 
numbers  are  gridpoints. 


GT.GID1  Grid  point  identification  numbers  (lnteger>0) 


Single-point  constraint  sets  must  be  present  in  the 
i nput. 

SPC  degrees  of  freedom  may  be  redundantly  specified 
as  permanent  constraints  on  the  GRID  card. 

A1 1  gridpoints  referenced  by  G1D1  thru  G1D2  must 
exi st. 

No  NASTRAN  "continuation"  cards  are  allowed. 


Field 

SID 

C 
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Input  Data  Card  ENDDATA 


Description:  Defines  the  end  of  the  Data  Deck 

Format  and  Example: 


1 

ENPDrtl  A 
ENDDATA 

2 

3  4  5 

6  7 

8  9 

10 

Remarks : 

1 . 

This  card  required 
exist  in  the  deck. 

even  if  no 

physical  data 

cards 

2. 

ENDDATA  nrjst  begin 

in  columns 

1  or  2. 

3. 

Failure  to  include  this  card 
operating  system  termination 
of  file  error. 

will  result  in 
caused  by  input 

an 

end 

1  * 


V. 


SAMPLE  PROBLEM 


This  sample  demonstration  problem,  taken  from  Reference  1 
is  presented  primarily  to  illustrate  the  new  NASTRAN  compatible 
data  input  procedure  and  format  included  in  OPTIM  III. 

The  structure  selected  is  a  symmetric  wing  box  as  shown  in 
Figure  4.  One  half  of  the  symmetric  ring  box  is  idealized  and 
consists  of  7  grid  points  and  16  finite  elements  of  the  following 
types : 

4  axial  force  members 
8  shear  web  elements 
2  quadrilateral  membrance  elements 
1  triangular  membrane  element 

Gridpoints  1  and  2  are  completely  constrained  in  the  x  and 
y  coordinate  directions.  Two  load  conditions  are  considered  and 
they  are: 

1.  5,000  Lb.  force  at  gridpoint  7  in  the  Z  coordinate 
direction. 

2.  10,000  Lb.  force  at  gridpoint  7  in  the  x  coordinate 
direction. 

A  summary  of  the  NASTRAN  compatible  input  data  (BULK  DATA) 
is  given  in  Figure  5.  Results  of  the  computations  performed  are 
given  in  Figure  6. 
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Figure  5  Input  Data,  Sample  Problem 
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L 
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74 
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26 
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27 
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25 

M 

20 

21 

52 

23 

24 
35 
34 
37 
51 
34 


PSFEAR 

PTRMEM 
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PMEE 

1 
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G 
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ft 
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I 
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a 
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Figure  6 


Output,  Sample  Problem 
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Figure  6  (Continued) 


INPUT  PROCESSOR  DIAGNOSTIC 
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BO  UN 

NO  ERRORS  OETECTED. 

_ OLOA 

_ NO.  ERRORS  DETECTED. _ 

ICON 

Md  ERRORS  DETECTED. 

NC  ERRORS  DETECTED 


ThE  COORDINATES  OF  THE  JiC POINTS 


O-C 

100.000 a 


70.0000 

7Q.COQO 


iniHHiiuti 


10.0000 


10.0000 

8.0000 


IllillUIllll 


lOO.OOOC 

iOO.OOOC 

140.C00O 

190.COOO 

8.0000 

3.0000 

Figure  6 

(Continued) 
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IN  ThE  X-CIRECT  ION 


IN  ThE  Y— C  IRECTION 


IN  ThE  Z-0 IRSCT  ION 


■***  WARNING  ***  SYMMETRIC  LGAQ  CONOITICN  I 

DEGREE  'GE  FREEDOM  I  IS  CCNSTRAI NED 

_ TUP  0.0 _ HAD  TS  SUPPR 

***  WARNING  ***  SYMMETRIC  LCAC  CCNDITICN  1 

DEGREE  CF  FREEDOM  2  IS  CCNSTRAI NED 


***  WARNING  ***  SYMMETRIC  LGAG  CQND ITICN  2 

.  DEGREE  OF  FREEDOM  I  IS  CONSTRAINED 

_ 1HJE _ LJ1 _ HAD  LS— SJi££Ri 

***  WARNING  ***  SYMMETRIC  LOAD  CONDITION  2 

DEGREE  CF  FREEDOM  2  IS  CCNSTRAI NEO 

_ IHE _ C*Q _ LCAC.  I.S S.UP.3R 

INCIV.CISPL  .CONSTRAIN  RED. DIR  NQS»  SYM  THEN  ANTI SYM 
ID 

- a - 

UPPER  ANC  LOWER  01 SP .CON STRA  1NT  LIMITS 
0.200  E  01 


0. 200  E  01 


0.0 

1CCOO.COCO 


Figure  6  (Continued) 
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Figure  6  (Continued) 


Figure  6  (Continued) 


. 


Figure  6  (Continued) 


V I .  PROGRAMMER'S  MANUAL 

The  programmer's  manual  for  the  OPTIM  III  Program  is  designed 
to  facilitate  its  implementation,  operation,  modification,  and 
extension.  This  manual  is  presented  in  the  appendix. 


VII.  CONCLUSIONS 


This  report  documents  the  revision  and  extension  of  the  weight 
optimization  program  OPTIM  II  to  a  new  program  OPTIM  III.  The  new 
program  utilizes  a  NASTRAN  input  type  format  (BULK  DATA  DARDS) 
and  is  applicable  to  the  analyses  and  optimization  of  complex  large  scale 
structures  which  contain  composite  membrane  plate  types  of  elements. 
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SECTION  I 


INTRODUCTION 


The  large  scale  automated  minimum  weight  design  program 
0PTIM  II  has  been  altered  to  make  its  input  format  compatible 
with  the  rlASTRAN  program.  The  resulting  new  program  is 
called  OPTIM  III  and  the  programming  aspects  are  described 
in  this  Appendix. 

The  information  presented  here  is  geared  to  the 
Programmer.  It  is  sufficient  to  fully  describe  the  program 
logic  and  the  required  peripheral  storage.  All  element 
generated  data  is  stored  externally  to  reduce  core  storage. 

A  separate  section  is  devoted  to  the  development  of  these 
files  so  that  I/O  time  may  be  optimized  through  efficient 
buffer  description.  Individual  subroutine  write-ups  are 
presented  along  with  the  complete  Fortran  source  listing. 

A  short  description  of  each  routine  is  included  to  aid 
in  getting  an  overall  familiarity  with  the  program  components. 
The  manner  of  allocating  storage  for  the  stress  and  displace¬ 
ment  matrices  is  described  in  diagram  form. 
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SECTION  2 

GENERAL  PROGRAM  LOGIC 

The  general  organization  of  the  "OPTIM"  program  is  illus¬ 
trated  in  Figures  1,  2  and  3. 


The 

program  consists  of  3  principal 

phases : 

1. 

Control  and  Data  Interpretation 

Phase 

2. 

Initialization  Phase 

3. 

Calculation  Phase 

1 .  Control  and  Data  Interpretation  Phase  (Figure  1) 

The  principal  purpose  of  this  pnase  is  to  read  the  NASTRAN 
form  input  and  prepare  the  information  for  the  next  two  phases. 
Five  subroutines,  MAIN,  SORT,  ZZ,  0PTIM2  and  OPIMPT  perf^^m  this 
function.  The  "MAIN"  -outine  controls  the  program  flow.  'he 
"SORT"  and  "ZZ"  routines  sort  and  count  the  data  cards  ana  then 
form  the  dynamic  sotrage  constants.  Each  NASTRAN  label  card  is 
processed  and  stopped  in  core  by  "ZZ".  An  input  file  data  tape 
based  on  OPTIM  related  input  is  then  written.  This  data  file  is 
then  used  by  "OPINPT"  which  writes  on  unit  NTAPE  all  the  data 
needed  for  initialization. 

2.  Initialization  Phase  (Figure  2) 

This  phase  produces  a  variety  of  information  for  the  Final 
Calculation  Phase.  Upon  being  activated  by  the  "MAIN"  routine, 
the  "NEWS'1  routine  determines  the  dynamic  storage  allocation 
for  the  Initialization  Phase.  The  "AONE"  routine  is  then 
activated  and  performs  the  following  tasks: 

A  load  matrix  is  generated  in  reduced  form  and  written  on 
tape  NSS 1 . 

Stiffness  and  stress  matrices  are  defined  by  the  called 
element  routines  {ELEM1  to  E L E M 8 )  and  are  written  cn  tapes  NSS1 
and  NSS2,  respectively.  The  element  routines  write  the  number 
of  stresses  in  the  element  and  the  weight  of  the  element  with  a 
unit  design  variable  on  the  tape  NSS3. 

A  final  tape,  NSS4,  is  written  with  information  on  linked 
groips,  stress  and  constraint  limits  and  various  other  key  control 
u  a  t  a . 

3.  Calculation  Phase  (Figure  3) 

The  Calculation  Phase  is  begun  by  the  "SIZE"  routine. 

Upon  being  calTed  by  the  "MAIN"  routine,  "SIZE"  determines  the 
'ynamic  storage  allocation  for  the  Incore  Calculation  Phase. 


For  the  case  where  there  is  enough  storage  available  for  an 
incore  solution  then  the  routine  "ATWO"  controls  the  calculations. 
If  enough  storage  is  not  available  then  control  is  passed  to 
the  "SIZE1"  routine  which  attempts  to  allocate  storage  for  an 
out-of-core  solution.  If  the  storage  for  out-of-core  is 
adequate  then  zhe  routine  "ATW01"  controls  the  calculation.  A 
brief  description  of  both  in-core  logic  and  out-of-core  logic 
follows : 

In  Core  Logic 

“The  routine  "ATWO"  controls  the  following  functions: 

A.  Tapes  MSS  3  and  NSS4  are  read.  These  tapes  are  read 
once  and  cne  information  is  used  to  initiate  the  program. 

B.  For  each  iteration  tapes  NSS1  and  NSS2  are  read.  The 
stiffness  matrices  are  assembled  in-core  and  then  displacements 
are  computed  by  the  banded  in-core  equation  solver  "BANCHO". 

All  the  element  stresses  are  computed  in  core  and  are  used  with 
the  displacements  by  the  ’’SCALE'*  routine.  Scaling  of  the 
structure  is  done  within  the  stress  and  displacement  constraint 
limits.  If  buckling  is  specified  then  during  the  scaling  process 
the  "BUCKTB"  and  "IGEN2"  routines  are  called  by  the  '.'SCALE* 
routine.  The  "MODIFY"  routine  then  adjusts  element  design 
variables  by  a  recursive  relationship  using  element  displacements 
and  stiffness  matrices. 

C.  The  print  out  of  displacements,  element  stresses  and 
element  design  variables  is  done  by  the  "ATWO"  routine. 

"ATWO"  also  calculates  and  prints  the  reactions  after  the  last 
iteration . 

Out-of-Core  Logic 

The  routine  "ATWOl,,controls  the  following  functions: 

A.  Tapes  NSS3  and  NSS4  are  read  and  used  to  initialize 
the  program.  Later  during  execution,  tape  NSS3  will  contain 
the  stress  matrix,  cne  column/record.  The  NSS4  tape  will 
contain  the  assembled  stiffness  matrix. 

B.  For  each  iteration  tapes  NSS1  and  NSS2  are  read.  The 
stiffness  matrices  are  assembled  out-of-core  on  tape  NSS4  by 
the  "BASS"  routine.  Tapes  NSS3  and  NSS5  are  used  as  scratch 
data  sets  for  this  step.  The  assembled  stiffness  matrix 
(which  is  also  banded)  is  triangularized  by  routine  "TCONTX" 
and  placed  on  tape  NSSJj.  Koutine  "EEQSX"  performs  the  final 
equation  solving  using  tape  NSS5  as  input.  The  stress  matrix 
is  computed  and  stored  on  NSS3  one  column/record.  The  stresses 
and  the  in-core  displacements  are  used  by  the  "SCALEl"  routine 
to  scale  the  structure  within  stress  and  displacement  constraint 
limits.  If  buckling  is  specified,  the  "SCALEl"  routine  calls 
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rwsas 


the  "BUCKTB"  and"lG£N2"  to  perform  this  function.  As  with  the 
in-core  solution,  the  "MODIFY"  routine  adjusts  element  design 
variables  by  a  recursive  relationship  using  element  displace¬ 
ments  and  element  stiffness  matrices  on  tape  NSS1. 

C.  The  print  cut  of  displacements,  element  stresses  and 
element  design  variables  is  done  by  the  "ATW01"  routine. 

"AOT01  also  calculates  and  prints  the  reactions  after  the  last 
iteration. 


Define:  Work  space 
logical  I/< 
control  va: 


Input  to 

Initialization 

Phase 


Utility  J 
Routines | 


w 


BCB 

MFGR 

GMPRO 

GTPRD 

MINV 


4 


-  Input  to 
Calculation 
Phase 


I  Element 
i  Matrix 
J  Generators 


2 (INITIALIZATION  PHASE) 
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SECTION  3 


EXTERNAL  DATA  SET  STRUCTURE 


This  program  uses  seven  data  sets  during  execution.  The  delivery 
version  of  the  OPTIM  program  comes  with  the  following  variable  names 
and  real  unit  designations.  These  unit  numbers  may  be  redefined  in  the 
"MAIN"  routine. 

Unit  Name  Unit  Id  Usage 

NTAPE/NSS5  8  Output  from  "MAIN"  and  "OPINPT" 

(BINARY)  routines  and  input  to  "AONE"  routine. 

For  out-of-core  solution  (Calculation 
Phase)  this  data  set  is  called  NSS5. 

It  is  next  used  as  a  scratch  data  set 
then  it  is  used  to  store  the  triangularized 
stiffness  matrix. 

NSS1  1  Contains  load  matrix  and  element 

(BINARY)  stiffness  matrices  with  assembly 

information 

NSS2  2  Contains  element  stress  matrices. 

(BINARY) 

3  ITOT  (no.  of  elements)  pairs  of  infor- 

( BINARY)  mation  from  element  matrix  generation 

routines  used  to  initialize  the  calculation 
phase. 

For  out-of-core  solution  (Calculation 
Phase)  this  data  set  is  next  used  as  a 
scratch  data  set  during  stiffness  matrix 
assembly.  Finally  it  is  used  to  store 
the  stress  matrix. 

4  Contains  initialization  information  for 

(BINARY)  ATWO  routine  or  ATWOl  Routine. 

For  out-of-core  solution  (Calculation 
Phase)  this  data  set  is  next  used  to  store 
the  banded  assembled  stiffness  matrix. 


-\  :  [ 
f 
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UNIT  NAME 


UNIT  ID 


USAGE 


L4 


9  Output  from  "SORT"  and  input 

(FORMATTED )  to  "ZZ".  Contains  identification 

record  and  NASTRAN  input  card 
image. 


L7 


7 

(FORMATTED) 


Output  from  "ZZ"  and  input  to 
"OPTIM2".  Contains  card 
images  in  OPTIM  MAGIC  format. 


NPIT/J5/L5  5  Standard  card  input  (SO  column  card) 

NPOT/  JE6/J6/L6  6  Standard  Line  Printer  (132  characters/ 

line). 


NOTE:  Since  units  5  and  6  are  standard  input  and  output,  they  will  not 
be  considered  in  this  discussion. 


NTAPE  CONTENTS 


OPDIVIR  Section  Input 

For  the  Out-of-Core  Calculation  Phase.  This  data  set  is  called  NSS5.  The 
triangularized  stiffness  matrix  is  stored  on  NSS5. 

NTCDS  20A4 

RECORDS 

Title  Cards 

One 

Record  ((Grid  (I,J),  J=1 ,  3),  1=1,  N2) 

(Coordinates) 

One 

Record 

Per  Element  NEL,  NID,  IBUCK,  (NODES(J),  J=1 ,  8),  (EM(J)J=1,  7) 

ITOT 

Records 

One  Record 

Linked  Group  IBGP,  (ILINK(LL) ,LL=1 ,IBGP) 

Bucking  Info 

ONE  RECORD 

Linkino  ( ILINK(LL) ,LL=1 ,  NSG) 

Info 

One  Record  (NBOUND(I),  1=1 , NAM)  NREACT  Array  number  of  each 

Boundary  component  bounded 

Condition 

Info 

NNZL 

Records  IR,IC,C1,C2 


Load 

Components 


IF 

NDL.NE.O 

One 

Record 

(NBDF(I),  1=1,  NDL) 

Individual 

Constraint 

Info  (ICON) 

IF 

NAA.NE.O 

NAA 

Records 

KNL,  (NBDF(JK) ,JK=1 ,KNL) 

Generalized 

Constraint 

Info  (GCON) 

If 

NAA+NDL.NE.O 

One 

Record 

(DISPU(I) ,DISPL(I) ,1=1 ,HDLL) 

IF 

IRST.EQ.l 

One 

Record 

(DVJR(J) ,J=1 ,IT0T) 

or  2 
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NSS1  Contents 


1st  Record 


3  Records 
Per  Element 


► 

(DELTA (I ),I=1,  Jl) 

i 

Load  Matrix 


J,K 


(X(L),L=1,J  ) 


Stiffness  Matrix 


(M(L) ,  L=1 ,  K) , (M2(L) ,L=1 , K) ,NNO , (NODES (L) , L=1 ,NNC } 


J  <- 


Symmetric  Antisymmetric 


Nodes  used  in 

Element 

Definition 


Assembly  info  vectors 
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NSS2  CONTENTS 


NSS3  Contents 


2  Records 
Per  Element 


LIST(I) 

^PAR(L)3L= 

KK,KK2) 

Number  of  Stresses 


Weight  Info  and 

info  for  eleme 


Purpose:  Initialization  of  A TWO  routine 


For  Out -Of -Core  Calculation  Phase 


The  stress  matrix  is  stored  on  this  data  set.  There  is  one 
record  for  each  column  of  the  stress  matrix. 
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NSS4  Contents 


NRDF , NLO AD , ITOT , NALD , NDL, NSG , NAA 


(SIGL(  I ) ,  SIGU ( I ) ,  ALL(  I ) ,  1=1,  ITOT) , 
IBS,  (iBK(I),  1=1,  IBS) 


Stress  and  Design  Paramet 
Limits  Buckling  Controls 


(N7(I)iI=l»NSG ) 


NRDF 


Linked  Group  Information 


(NBDF(I),I=1,K),(NBDF2(I),I=1,K) 


No.  of  Reduced  Deg. 
of  Freedom 


Individual  Constraint  Inf 


NEAA,  ( NE ( I ) , 1=1 , NEAA} ,  NE  A2 &1Z2 (I),  1=1, NE A2 ) 


Generalized 
Constraint  Info 


( DISPU ( I ) , DISFL(I) , 1=1,  NDL) 


Limits  for  Constraints 


Cl, C2 , ITERN 


( AREA( I ) ,  1=1, ITOT 


L0AD1,  L0AD2 ,  NSYM ,  N  ASYM,  NONSTM,  NRDF, 

( LNOD( I ) , 1=1 , NRDF ) , NRDF2, ( LN0D2 (I) ,  1=  I,  NEDF2 ) 


NREF, NDOF,KKTROL, NBOU, ( NPI ( I ) , 1=1, NBOU) , 
NB0UN2 , (NP2(I) ,1=1, NB0UR2 ) 


For  Out-Of-Core  Calculation  Phase 


The  banded  and  assembled  stiffness  matrix  is  stored  in 
this  data  set. 


L4  Contents 


(80  character  card 
images  > 


These  records  can  be  in  any  order  on  the  file.  It  represents  data  in  "pairs" 
of  records.  Each  pair  represents  one  specific  type  of  data  items  determined 
by  the  value  "NKIND"  and  "NCOUNT"  which  are  written  on  the  first  record 
of  the  pair.  R 

Record  one  of  the  pair:  NKIND,  NCOUNT 

Record  two  is  determined  by  NKIND  as  shown  in  the  following  table: 


NKIND 


Second  Record 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


BULK  data 

TITLE  " 

BUCK 

CONROD  " 

COD  MEM 

CQDMEM1 

CQDMEM2 

COD  MID 

CQUAD1 

CQUAD2 

CRD  MID 

CROD 

CSHEAR 

CTRIA1 

CTRIA2 

CTRMEM 

CTRMID  " 

CTUBE 

END*  " 

ENDDATA 

FORCE  " 

GCON 

GRID 

ICON 


* 


25 

LINKS 

data 

26 

MOMENT 

tt 

27 

MAT1 

ti 

28 

OPLOAD 

tt 

30 

OPDVIR 

i 

31 

OPTIM 

i . 

32 

PQDMEM 

it 

33  ■ 

PQDMEM1 

it 

34 

PQDMEM2 

it 

35 

PQDMID 

n 

36 

PQUAD1 

it 

37 

PQUAD2 

tt 

38 

PRD  MID 

ti 

39 

PROD 

tt 

40 

PS  HEAR 

ii 

41 

PTRIA1 

it 

42 

PTRIA2 

ii 

43 

PTRMEM 

ti 

44 

PTRMID 

it 

45 

PTUBE 

it 

46 

SPC 

it 

47 

SPC1 

tt 

48 

SPC1  THRU 

tt 

49 

CWEB 

it 

^0 

PWEB 

tt 

51 

BEGIN  BULK 

tt 

52 

CTRIM6 

it 

53 

CTRSHL 

tt 

54 

CQUAM8 

it 

55 

CQUSHL 

tt 

56 

PTRIM6 

it 

57 

PTRSHL 

it 

58 

PQUAM8 

it 

59 

PQUSHL 

it 
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L7  Contents 


(80  character 
card  images) 


Record  #1 
"  #2 
"  #3 

"  #4 

"  #5 

"  #6 
"  #  7 

"  #8 


NGRI 

records 

NEL 

records 


NECARD 

records 


NOCARD  records 


N LOADS 
TIMES 


NGR 

records 


NICON 

records 


REPORT  (2A4) 

TITLE  (2A4) 

NTIT  (6X,  113) 

(IEDING(I),  1=1,  20) 

SYSTEM,  (BLANK,  1=1,  18) 

NGR1,  NGR2,  NS4S,  ILOAD,  NECARD,  MNN 
OPTIM  (OPT (I),  1=1,  16) 

COORD,  (BLANK,  1=1,  18) 

■£m.  Cl,  C2,  C3 
ELEM,  (BLANK,  1=1,  18) 

£l,  NOPT(I),  X,  NNODES(I),  (NODENO(I,  J),  J-l,  8) 
OEXTERN,  (BLANK,  I+i:  18) 

fl,  BUCK(I),  A,B,C,D,E,F,G,H,P,Q,R, 

T  J,  ANGLE(I) 


LINK,  (BLANK,  1=1,  18) 
LLINK 

BOUND,  (BLANK,  1=1,  18) 
(BLANK,  1=1,  20) 

MODAL  (1,  1=1,  NSIS) 

I,  (LBOUND  (I,  J),  J  =  l,  NSYS) 

OLOADS,  (BLANK,  1=1,  18) 
LCOND, IL 

MODAL,  0,0,0,0,0,0,0,40,0 
IGR,  Cl,  C2 ,  C3,  C4 


ICON,  (BLANK,  1=1,  19) 
^EYEC  Cj.  I).  J  =  l.  13) 


N GROUP 
TIMES 


NGROUP 
J  records 


{ 


GCON,  (BLANK,  1=1,  19) 
GGROUP,  I 

GROUP (5,  I),  GROUP (6,  I) 


NGC1 
\  record 


{GCOND  (?•,  J),  (GCOND(K,  J)  K=3,  14) 


.  OPDVIR,  (BLANK,  1=1,  18) 


NLL 

fOPD,  (BLANK,  1=1,  18) 

TIMES 

\0PDV 

last 

FEND  =  "END" 

record 

1 


Section  4 

Subroutine  Write-Ups 


The  OPTIM  program  contains  37  subprograms ,  each  with  a 
unique  sequence  number.  Columns  73,  74,  75  and  76  contain  the 
"DECK"  name  and  columns  77  through  80  contain  the  card  sequence 

number  for  that  subprogram.  The  first  is  always  - 0000  with 

successive  increments  of  10. 

Included  with  each  description  is  a  statement  declaring 
the  size  of  the  subprogram.  This  number  is  intended  as  a  guide 
only,  as  it  reflects  the  storage  requirement  on  an  IBM/360/65, 
FORTRAN  G  level  19  compiler. 

In  this  manual  the  subroutine  writeups  are  presented  in 


alphabetical  order. 


Brief  Routine  Descriptions 


Routine 

Name 

ADJUST 

AONE 

ATWO 

ATWOl 

BANCHO 

BASS 

BOB 

BEQSX 

BMAT 

EUCKTB 

CMAT3 

CUBIC 

EGMAT 

EMAT 

ELEM1 

ELEM2 

ELEM3 

ELEM4 

ELEM5 

ELEM6 

ELEM7 

ELEM8 


Purpose 

Adjust  each  NASTRAN  form  field  so  that  it  is  either 
left  or  right  adjusted. 

Control  initialization 

Control  calculation  and  print  out  (in-core  solution) 
Control  calculation  and  print  out  (out-of-core  solution) 
Banded  in-core  equation  solver 

Assemble  banded  stiffness  matrix  (out-of-core  solution) 
Matrix  triple  product 

Equation  solver  (out-of-core  solution) 

Generates  A  matrix  for  triangular  plate,  also  area 

Interfaces  with  SCALE  and  IGEN2.  Determines  minimum 
plate  thickness  for  stability. 

Generates  C3  matrix  for  triangular  plate 

Solve  for  principle  beam  stresses 

Generates  EG  matrix  for  triangular  plate 

Generate  elasticity  matrix  for  either  orthotropic  or 
isotropic  properties ,  membrane  triangle 


'S 


IElemert 
Matrix 
Generators 


Axial 
Shear  web 
Triangle 
4  Quad 

Beam  (tube) 
Mid-point  axial 
Mid-point  triangle 
Mid-point  quad 


Routine 

Name 

ESCONT 


FOMO 


JMPRD 


GTPRD 


IGEN2 


INSPC 


INSRTI 


KCALC 


KMAT 


MAIN 


Brief  Routine  Descriptions  (Cont'd.) 


Purpose 

Solve  matrix  equation  A*  A  Transpose  *  X  =  F,  for  X 
(out-of-core  solution) 

Process  Force  Cards 

Matrix  product 

Matrix  product  with  transpose  of  first  matrix 

Solves  Eigenvalue  for  buckling  coefficient  of  plate 

Stores  data  obtained  from  input  card  into  working 
storage 

Stores  element  properties  and  material  numbers  into 
working  storage 

Solve  for  XK  where  A  *  XF  =  F  and  A  is  banded  lower 
triangular  matrix  (out-of-core  solution) 

Generates  stiffness  matrix  for  triangular  membrane 

Control  program  flow 


a«  f 

P  ! 


MFGR 

MINV 

MODIFY 

NEWS 
OPT I M2 
OPINPT 


PS3R 


READ  I 


SCALE 


SCALE 1 


SIZE 


Matrix  Ranking 


Matrix  Inverse 


Next  element  design  parameter  guess  after  scaling 

Allocate  storage  for  initialization 
Interprets  input  and  initializes  data. 

Interpret  report  from  input  and  place  results  on  NTAPE 

Define  cubic  equation  from  beam  stress  components 

Normalize  stresses 

Reads  and  modifies  input  data 

Scale  structure  using  stress,  displacement  and 
buckling  constraints  (in-core  solution) 

Scale  structure  using  stress,  displacement  and 
buckling  constraints  (out-of-core  solution) 

Allocate  storage  for  calculation  phase  (in-core 
solution) 


Brief  Routine  Descriptions  (Cont’d) 


Routine 

Name  Purpose 

SIZE1  Allocate  storage  for  calculation  phase  (out-of-core 

solution) 

SMA2  Generate  stress  element  matrix  for  membrane  triangular 

plate  with  orthotropic  properties 

SORT  Sort  and  count  data  based  on  LABEL  information 

SPCSOB  Process  SPC  (single  point  constraint)  cards 

TCONTX  Control  tape  flow  to  triangularize  stiffness  matrix 

(out-of-core  solution) 

TilAT  Generate  transformation  matrix  for  triangular  plate 

orthotropic  angle 

TTRI  Triangularize  rows  of  a  banded  matrix  A  (out-of-core 

solution) 

WRITEL  Tests  character  of  element  connection  card  and  writes 

element  information  on  to  file 


i.4 


XCALK 

XTRAK 


Solve  for  x  where  A  *  X  =  XK  and  A  is  an  upper 
triangular  matrix 

Interpret  degree  of  freedom  informations 


1 

a 


XTRAK2  Interpret  degree  of  freedom  informations  for 

generalized  constraint  information 

ZZ  Generates  OPTIM  data  which  is  input  by  NASTRAN  form 

input  cards 


jf 

aM 


i 


Subroutine  Name 


ADJUST 


Purpose 

Procedures 

Input  Arguments 


Output  Arguments 

Error  Returns 
Calling  Sequence 
Subroutine  User 


Adjust  each  NASTRAN  form  field  so  idiat 
it  is  either  left  or  right  adjusted. 

Each  character  of  the  word  is  tested  to 
see  if  it  is  blank.  When  a  non-blank 
is  met,  it  is  shifted  to  the  end  of  the 
word. 

WORD  -  Input  word  to  be  either  right 
or  left  adjusted 

Icode  Icode  =  8  specifies  that  the 

word  should  be  right 
adjusted 

Icode  =  1  specifies  that  the  word 
should  be  left  adjusted 

Word  right  adjusted,  word  is  stored  back 
into  WORD 

None 

CALL  ADJUST  (Word,  ICODE) 

NEWS 


Subroutine  Name: 


AONE 


Purpose:  Assemble  all  initialization  information 

for  the  calculation  phase  of  the  program. 


Equations  and  Procedures : 

a)  Print  out  title  cards. 

b)  Retrieve  (15)  coordinate+element  data.  If 

requested  print  this  data  out. 

c)  Recalculate  boundary  conditions  if  symmetry 

plane  nodes  are  specified. 
d|  Create  load  matrix  from  input  loads. 

e)  Define  virtual  loads  for  individual  constraints 

and  generalized  constraints. 

f)  Call  element  routines. 

g)  Pour  tapes  of  data  are  written  with  information 

for  the  calculation  program  (ATWO). 

Input  Arguments : 


15 

SPOT 

SSS1-NSS4 

N2 

SAA 

NDL 

NSG 

NNZL 

KLN2 

NSYM 

NASYM 

NO  NSYM 

ITERN 

NTCDS 


(Unit  8)  tape 
(Unit  6;  printer 
(Units '  1-4 )  tapes 
No.  of  nodes 

No.  of  generalized  constraints 

No.  of  Individual  Constraints 

No.  of  Linked  Groups 

No.  of  input  loads 

No.  of  symmetry  plane  nodes 

No.  of  summetric  load  conditions 

No.  of  antisymmetric  load  conditions 

No.  of  non-symmetric  load  conditions 

Max.  no.  .of  iterations 

No.  of  title  cards 


Output  Arguments : 

Convergence  limits 

No.  of  reduced  DOF  (symmetric  load  condition) 
No.  of  generalized  displacement  constraints 
Total  no.  of  constrained  DQF  (symmetric 
load  condition) 

Calculation  control 

Error  Returns :  IER 

Calling  Sequence: 

CALL  AONE (X, Y,Z,N5,N6,N7,E,AMU,N8,N17,ALL, RHO,SIGU, 

S IGL. R,  Nil, N13,N15, NOAL, N0AL2, KL, KL2,  LNOD,  LNOD2, 
DISPUjDISPL,  NBDF, NBDF2, NSE, DELTA ,15,  KPOT, NSS 1, NSS  2 , 

NSS 3,N3S4, IBUKL,  IELT,  I3K, Nl, N2,NAA. NDL, NSG, NNZL, 


C1INP) 

C2INP) 

NRDF 

NEAA 

NBOU 

IRST 


NALD,  XLN2,NSYM, NASYM,  NONSYM,  ITERN,  NTCDS .  E5SS1 ,  ITOT, 
IELI,  NREACT,  C IMP ,  C  2IHP ,  NRDF,  NRDF2 ,  NEAA ,  NEA2 ,  NDL1 , 
NDL2,  NBOU,  NB0UN2,  IRST,  SNCAL,  ISR ) 

8.  Input  Tapes.:  NTAPE  (unit  8) 

9.  Output  Tapes: 

NSS1  (unit  I) 

NSS2  (unit  2) 

NSS3  {unit  3) 

NSS4  (unit  4) 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (20160  bytes)  5040  words 

12.  Subroutine  User:  NEWS 

13.  Subroutine  Required: 

ELEM1 

ELSM2 

ELEM3 

ELEM4 

ELEM5 

EEEM6 

ELEM7 

SI2M8 

14 .  Remarks :  None 
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1.  Subroutine  Name: 


A  WO 


2.  Purpose:  Control  execution  of  the  in-core  calculation 

phase  of  the  program. 

3.  Equations  and  Procedures: 

a.  Read  in  all  data  from  NSS4  tape. 

b.  Generate  information  vectors  for  BANCEO  (in-core 

equation  solver). 

c.  Read  in  data  from  NSS3.  Element  information 

consisting  of  number  of  stresses  and  weight 
factor. 

d.  Assemble  stiffness  matrices  and  solve  for  stresses 

and  displacements  in  two  passes.  The  first  pass 
computes  symmetric  load  conditions  and  the 
second  pass  computes  anti -symme trie  load 
conditions . 

e .  If  this  is  a  statics  problem  go  to  g)  otherwise 

continue . 

f.  Call  SCALE  routine  to  scale  stresses,  displacements 

and  design  variables . 

g.  Print  design  parameters,  displacements  and  stresses 

if  requested  for  every  iteration  or  if  this  is 
last  iteration.  If  statics  run  go  to  j). 

h.  Check  convergence. 

i.  Call  MODIFY  routine  and  return  to  d). 

j .  If  last  iteration  compute  and  print  reactions  by 

load  and  node. 

4.  Input  Arguments: 

LLX  Total  number  of  stresses 

NPCT  (Unit  6)  printer 

NSS1-NSS4  (Units  I~ii)  tapes 

NGV  Calculation  control 

NT3RN  Maximum  number  of  iterations 

IPi  Print  control  for  displacements 

I?2  Print  control  for  stresses 

IPRINT  Print  control  for  element  design  variables 

I3UEL  Buckling  control  from  OPTIM  section  of  input 

5.  Output  Arguments:  None 

6.  Error  Returns :  IE?. 

7.  Calling  Sequence: 

CALL  A  TWO  (NT,  SIGL,  SIGU,  ALL,  LIS  T,  AREA,  AQEA,  3IGS ,  PAR,  DELTA, 
STRESS ,  NE ,  NE2 ,  DISRJ,  DISPL,  If  PDF,  NBDF2 ,  LNOD ,  LA ,  LR.  AUX, 

N2  ,LA  2 ,  LNCD2,  ST ,  NPI ,  N?2,ADUA ,  E  5  ,  N6 ,  AREAD ,  LLK,  SPOT,  NSS 1, 
NSS2,NSS3,NSS4,  IPI,  I?2,NGV.  STERN,  IPRINT,  LBUXL,  I3K,  IER) 
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3 .  Input  Tapes : 

NSS1  (Unit  1) 

NSS2  (Unit  2) 

NSS3  (unit  3) 

NSS4  (Unit  4) 

9.  Output  Tapes:  None 

10 .  Scratch  Tapes :  None 

11.  Storage  Required: 

(14032  Bytes }  3508  words 

12.  Subroutine  User:  Size 

13.  Subroutine  Required: 

BANC HO 

SCAI2 

MODIFY 

14 .  Remarks :  None 
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1.  Subroutine  Name:  ATW01 


2.  Purpose:  Control  execution  of  the  out-cf-core 

calculation  phase  of  the  program. 

3.  Equations  and  Procedures: 

The  difference  between  the  AIWO  and  A3W01  routines  is 
that  the  stress  matrix  is  stored  on  tape  and  the 
assembled  stiffness  matrix  is  also  on  tape  for  the  ATW01 
routine . 


a)  Read  in  all  data  from  NSS4  tape. 

b)  Read  data  from  NSS3  tape. 

c)  Assemble  stiffness  matrices  and  solve  for  displace¬ 

ments  and  stresses  in  two  passes.  The  first  pass 
selves  for  symmetric  load  conditions  and  the 
second  oass  solves  for  anti-symmetric  load 
conditions.  Call  routines  BASS,  BEQSX  AND  TCONTX 
to  assemble  and  solve  for  displacements. 

d)  Go  to  f)  for  statics  problem. 

e)  Call  SCALE1  routine  to  scale  stresses,  displace¬ 

ments  and  design  variables . 

f)  Print  design  variables,  displacements  and  stresses 

if  requested  for  every  iteration  or  if  this  is 
the  last  iteration.  If  statics  run  go  to  i). 

g)  Check  con verge  nee  . 

h)  Call  MODIFY  routine  and  return  to  c ) , 

i)  If  last  iteration  compute  and  print  reactions 

by  load  and  node . 


4.  Input  Arguments: 


LLK 

NPOT 

NSS 1 -NSS 5 
NGV 
NTERN 
IP1,I?2,  ) 
IPRLNT  ) 
IBUKL 


Total  number  of  stresses 
Printer 
Tape  Units 
Calculation  control 
Maximum  number  of  iterations 
Print  controls  for  displacement  stress+ 
element  design  variables,  respectively 
Bucxling  control  from  OPTIM  section  of 
input 


5.  Output  Arguments:  None 

6.  Error  Returns:  None 


7.  Calling  Sequence: 

ATW01  (N?,SIGL,SIGU,  ALL,  LIST,  AREA,  AQEA,  BIGS,  PAR, 
DELTA , 3 'TRESS ,  NE ,  NE2,DI3?U,DISPL,  NBDF,  NBDF2,  LA,NZ,LNOD2, 3T, 
NP1,  NF2 , ADEA , N5 . N6 , AREAD , LLK, NPOT, NSS 1, NSS  2 , NSS 3 , 
N3S4,NSS5, IP1, IP2,NGV, NTERN, IPRINT, IBUKL, IBK, NTOTAL, 

ISR)  . 


I 
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8 .  Input  Tapes : 

NSS1 

NSS2 

NSS3 

NSS4 

NSS5 

9.  Output  Tapes:  Non 

10.  Scratch  Tapes:  NS-C 

11.  Storage  Required: 

(14978  Bytes)  3745  voids 

12.  Subroutine  User:  S'JZ 

13.  Subroutine  Required: 

BASS 
SCALE  1 
MODIFY 
BEQSX 
TCONTX 


1. 


Subroutine  Name : 


BANC  HO 


2.  Purpose: 

In  core  banded  matrix  equation  solver. 

3.  Equations  and  Procedures: 

When  given  an  N by N positive  definite  matrix  [a]  a 
general  CholesXy  triangularization  is  effected. 

4.  Input  Arguments: 

A  Matrix  of  coefficients 

R  Independent  variable  (vectors) 

N  Number  of  equations 

M  Number  of  dependent  variables  for  solution 

IR  Row  numbers 

IC  Column  number  of  first  non-zero  term  in  A  (for 
each  row) . 

IA  Diagonal  locations  in  A 

AUX  Dummy  array  (N  long) 

5.  Output  Arguments: 

R  Solution  vectors 

A  Is  destroyed 

6.  Error  Returns: 

Writes  matrix  can  not  be  inverted,  ROW,  l4,  STOP 

7.  Calling  Sequence : 

CALL  BANC HO  (A, R, N,M, IR, IC, IA, AUX ) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required: 

(4358  Bytes )  IO89  words 

12.  Subroutine  User:  ATWO 

13.  Subroutine  Required:  None 
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1. 


Subroutine  Name : 


BASS 


i 


t 

I 


i- 


I 


■j 
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2.  Purpose; 

Assembly  of  banded  symmetric  matrix  from  symmetric 
element  arrays. 

3.  Equations  and  Procedures: 

aj  The  total  assembled  matrix  is  initially  set  =  zero. 

b)  A  "core  full"  section  of  the  assembled  matrix  is 

handled  at  one  time.  Each  element  matrix  is  read 
into  STA  from  tape  14.  Each  LISTEL  is  read  from 
tape  14.-  This  LISTEL  is  an  element  matrix  array 
which  contains  the  row  numbers  of  the  assembled 
matrix  into  which  the  element  matrix  must  be 
added . 

c)  Each  fixed  point  number  in  LISTEL  is  examined.  If 

it  Is  =  0  no  addition  into  the  system  matrix  is 
made.  If  it  is  greater  than  the  maximum  row  number 
considered  in  the  current  computation  pass,  the 
particular  element  information  is  put  on  Tape 
18  for  later  use  in  a  subsequent  computation 
pass . 

d)  If  more  than  one  computation  pass  is  necessary, 

the  section  of  the  element  matrix  which  has  not 
been  used  is  placed  on  tape  13. 

e)  Steps  2,  3,  4,  are  repeated  for  each  element. 

f)  The  assembled  matrix  portion  for  each  computation 

pass  is  placed  on  tape  13  in  banded  form.  Each 
row  uf  the  matrix  constitutes  a  record. 


4.  Input  Arguments: 

STA  Storage  array  for  assembled  matrix 

NELEM  Number  of  elements 

NMDB  Order  of  system 

LISTEL  Element  of  assembly  vector 

NZEL  Cumulative  total  of  non  zero  elements  in 

rows  1  through  i-1  or  reduced  matrix 
IZr  Number  of  zero  elements  In  row  of  reduced 

matrix 

IRA  Work  storage  array 

IMAX  Number  of  work  storages  available 

ST  Storage  array  for  element  stiffness  matrix 

AREA  Element  design  parameter  array 

KASEM  Symmetric/anti-symmetric  indicator 
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Output  Arguments : 

N5  Output  element  array  read  from  tape  l4 

N6  Output  element  array  read  from  tape  14 

BIGS  Output  element  array  read  from  tape  14 

Error  Returns:  None 

Calling  Sequence: 

CALL  BASS  (STA,NEIEM,  NMDB,LISTEL,  NZEL,  IZR,  IRA,  IMAX, 
ST,l4,l8,I3,I7,ArJEA,KASEM,N5,N6,BIGS). 

Input  Tapes : 

Tape  14  -  Element  matrices  and  element  controls 
Output  Tapes : 

Tape  13  -  Assembled  matrix  in  banded  symmetric  form. 
Each  row  equals  a  record. 

Scratch  Tapes : 

17  -  Temporary  storage  for  element  matrices 

18  -  First  temporary  storage  for  element  matrices 

for  later  computation  pass 

Storage  Required: 

(3456  Bytes)  864  words 
Subroutine  User:  ATW01 


Subroutine  Required:  None 


Subroutine  Name: 


BCB 


Purpose:  To  evaluate  the  triple  product:  the . transpose 
of  a  matrix  A,  a  symmetric  matrix  S  and  the  A  matrix 

Equavlons  and  Procedures: 

T 

ANMM  =  H  £1  AMN  *  SNN  *  anm 
N  N 

Input  Arguments : 

A  The  elements  of  the  A  matrix 

SYM  The  elements  of  the  S  matrix  (symmetric  - 

bottom  half) 

ND,MD  Dimension  of  A  matrix 

N,  M  Order  of  A  matrix 

N1  Number  of  rows  to  be  deleted  in  multiplication 

=  0 

SCAL  Scalar  quantity 

IASSY  See  remark 

Output  Arguments : 

AN:  Elements  of  the  matrix  AN  which  is  the  final 

product 

Error  Returns:  None 

Calling  Sequence : 

CALL  ECB(A,  SYtf, AN, ND,  MD , N,  M,  Nl,  SCAL,  IASSY ) 

Input  Tapes :  None 

Output  Tapes :  None 

Scratch  Tapes :  None 

Storage  Requir  ed : 

(1040  bytes)  26u  words 
Subroutine  User: 

ELEM7,  ELEM4,  ELEM5 
Subroutine  Required:  None 
Remarks : 

IASSY  controls  the  summation  procedure. 

If  =  1,  AN  will  be  the  sum  of  the  calculated  AN  and  all 

previous  elements  of  AN. 

=0,  AN  will  be  the  triple  product. 


1. 


Subroutine  Name : 


beqsx 


2.  Purpose: 

To  perform  Simula tanecus  aquation  solution  for  banded 
symmetric  matrix  input  using  Cholesky  procedure. 

3.  Equations  and  Procedures: 

1.  The  program  is  designed  so  that  a  "core-filled" 

piece  is  considered  at  one  time.  The  procedure 
to  handle  this  is  set  up. 

2.  All  loads  are  stored  in  pool  from  the  input  argument 

list. 

3.  A  call  to  ESCONT  is  made.  This  is  the  routine  which 

actually  does  the  computation. 

4.  Displacements  are  stored  in  the  DISPL  array  output 

from  ESCONT. 

4.  Input  Arguments: 

AK  Storage  array  for  input  matrix 

NZEL  Cumulative  total  of  non-zero  elements  in 

rows  1  through  i-1  of  reduced  matrix 
IZR  Number  of  zero  elements  in  row  of  reduced 

matrix 

XK  Working  storage  array 

PCOL  Load  vector 

NMDB  Order  of  system 

NTOTAL  No.  of  work  storages  available 

NL  No.  of  load  conditions 

5 .  Output  Arguments : 

DISPL  -  Displacements 

6.  Error  Returns:  None 

7 .  Ca lling  Sequence: 

CALL  BEQSX  (AK,NZEL,  IZR, XK, DISPL, PCOL, I7,NMDB,  NTOTAL, NL) 

8.  Input  Tapes: 

17  -  Input  triangularis ed  banded  matrix  array 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required: 

(1864  Bytes)  466  words 
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Subroutine  NAME: 


2.  Purpose 


3 .  Equation  and 
Procedures 


4 .  Input  Arguments 


5 .  Output  Arguments 


6 .  Error  Returns 


7 .  Calling  Sequence 


8.  Subroutine  User 


BMAT 


Generates  A  matrix  for  triangular 
plate.  Also  generates  AREA  of  triangle. 


BINV  _ 


x2y3 

y2 


x2y3 


x3y2 

y3 


x3y2 


y3  -*2 


Area  —  x2^3  x3^1  —  ^r3xl  —  x3^2 


x  (i) 
y  (i) 

BINV 

AREA 

None 


x2yl> 

x  coordinate  values  1  ^ 
y  coordinate  valuesj 


=  1,  2,  3 


output  matrix 
output  value  of  area 


CALL  BMAT  (BINV,  X,  Y,  AREA) 


ORTHO 3 


9 .  Subroutine  Required  None 


Subroutine  Name : 


3UCXT3 


f 


I 


i 


Purp  ose : 

This  subrcuti .e  will  determine  the  minimum  plate 
thickness  for  the  stabiliry  of  a  flat  rectangular 
simply  supported  plate  subjected  to  biaxial  stresses 
and  shear. 


3- 


Equations  and  Procedures : 


The  Shear  3uckling  Coefficient,  Xg,  is  determined  for 
symmetric  and  antisymmetric  buckling  modes  wherein  the 
eigenvalue  subroutine  IGEN2  is  employed.  A  one 
dimensional  Newton  F.aphson  iterative  procedure  is  used 
to  converge  the  K  calculated  and  the  K  determined 
from  the  element  geometry  and  state  of  stress.  The 
solution  is  based  on  the  Raleigh -Hi  tz  method  and  the 
general  form  assumed  for  the  deflected  surface  is 


A  sin 
mn 


oXi"  x 
a 


sin 


nTTj- 

b 


Xs  is  calculated  from  the  homogeneous  linear  set  of 
equations  represented  as 


mnna 

(m2-?2) (n2-q2 ) 


ana 


■rf  2 

52^7? 


4. 


where  fi=  a/b;  m=l,  2,  3  i  n  =  1,2,3  • and 

(m±p,n±q)  are  cdd  numbers  . 


Input  Arguments : 


3IGX 

SIGY 

TAU 

T0 

A L,  3L 
DC 

ISLE 


Stress  in  the  x  direction 
Stress  in  the  y  direction 
Shear  stress  in  x-y  plane 
Given  thickness 

Lengths  of  sides  of  rectangle  AL  31- 
Input  material  constant 
Element  number 


Output  Arguments; 

T3  -  Output  material  thickness 


1  U 


Calling  Sequence: 

CALL  3UCXTB ( S IGX, SIGY, TAU , TO . AL, SL,DC , ISLE, IOP , TB ) 

Input  Tapes :  None 

Output  Tapes :  None 

Scratch  Tapes :  None 

Storage  Required: 

(9880  Bytes)  2^70  words 
Subroutine  User:  SCALE 

Subroutine  Required:  IGEN2 


1 .  Subroutine  NAME 


CMAT3 


I  2.  Purpose 

I 

3.  Equations  and 

Procedures 

: 

i 

t 

( 

t 

i 

I 

* 

4 .  Input  Arguments 

5.  Output  Arguments 

6 .  Error  Return : 

7.  Calling  Sequence 

8.  Subroutine  User 


Generates  matrix  for  triangular  plate 
=  T  *  Eg  *  TT 


T  = 


0 

0 

0 

0 

0 

1 


EG 


3x3  material  elasticity  matrix 


3 

None 


Matrix  defined  by  equation  above 


CALL  CMAT3  (C7  EG) 
J1 


ORTHO 3 


1.  Subroutine  Name:  CUBIC 

2.  Purpose: 

To  determine  principle  stresses  for  the  tube  (beam) 
element . 

3.  Equations  and  Procedures: 

Where  X  is  the  stress  given 
X3  +  AX2+BX+C  =  0 

Solve  the  cubic  equation  and  return  the  three  roots. 

4 .  Input  Arguments : 

A,  B ,  C  -  Coefficients  of  cubic  equation 
XX  -  Proportionality  factor 

5.  Output  Arguments: 

A,B,C  -  Three  unequal  real  roots 

6.  Error  Returns:  None 

7.  Calling  Sequence:  GALL  CUBIC  (A,B,C,XX) 

8 .  Input  Tapes :  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required: 

(1748  Bytes)  437  words 

12.  Subroutine  User:  PS3D 

13.  Subroutine  Required:  None 

14 .  Remarks :  None 


1 .  Subroutine  Name 


EGMAT 


2.  Purpose  Generate  Eg  matrix  for  tx.'  angular 

plate 

3.  Equation  and  Eg  =  TET1 

Procedures 

4.  Input  Argument  E  =3x3  matrix 

T  =3x3  transformation  matrix 

5.  Out  Arguments  Eg  =  matrix  defined  by  equations 

above 

6.  Error  Returns  None 

7.  Calling  Sequence  Call  EGMAT  (EG.  E,  T) 

8.  Subroutine  User  0RTH03 
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1.  Subroutine  Name 


EMAT 


2.  Purpose  To  generate  elasticity  matrix  for  either 

orthotropic  or  isotropic  properties 


3.  Equations  and 
Procedures 

Orthotropic 
(IMAT  =  2) 

E 

= 

G11 

G12 

G13 

G21 

G22 

G23 

G31 

G32 

G33 

Isotropic 
(XMAT  =  1) 

E 

=| 

Ex 

Muxy 

0 

MxyEy 

Ey 

0 

1 

0 

0 

Gxy 

4.  Input  Arguments  Ex,  Ey,  M  ,  G  -  material  properties 

xy  xy 

EM  -  contains  values 

IMAT  -  Orthotropy  code  =  2  orthtropic 

Orthotropy  code  =  1  isotropic 

5.  Output  Arguments  E  >-  output  3x3  matrix  defined  by 

equations  above 

6.  Error  Returns  None 

7.  Calling  Sequence  Call  EMAT  (E^,  Ey,  M^,  G^,  E,  EM,  IMAT) 

8.  Subroutine  User 
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1 .  Subroutine  Name :  ELEM1 


2.  Purpose: 

To  compute  the  elemental  stress  and  stiffness  matrices 
for  an  axial  element. 


3 .  Input  Arguments : 


X,  Y,  Z 
E 

AMU 

RHO 

N5,  N6 
NEL 
N1 
N2 


Node  point  coordinates 
Modulus  of  elasticity 
Poisson's  ratio 
Density 
Element  nodes 
Element  number 
Degrees  of  freedom/node 
Total  number  of  nodes 


4.  Output  Arguments: 


M 

M(25) 

M(26)=1 

ST(l') 

ST(301) 

AK(1) 


Unreduced  direction  numbers  for  degrees  of 
freedom  of  element 
Order  of  stiffness  matrix 
Number  of  stress  components  (tape  NSS3) 
Start  of  stiffness  matrix 
Start  of  stress  matrix 

(on  tape  NSS3)  weight  of  element  with  unit 
design  variable 


5.  Error  Returns:  None 

6.  Calling  Sequence: 

CALL  EIEM1  (X,  Y,Z  E,AMU,RH0,N7,N5,N6,N8,N9,N10,N11, 
N12,R,M,NELj  N1,N2) 


7. 

Input  Tapes: 

None 

8.' 

Output  Tapes : 

NSS3 

9- 

Scretdh  Tapes: 

None 

10. 

Storage  Required: 

(1789  Bytes)  445  words 

11.  Subroutine  User:  AONE 

12.  Subroutine  Required:  None 

13 .  Remarks : 

The  stress  matrix  provides  one  stress  component, 0"x. 
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1.  Subroutine  Name: 


ELEM2 


2.  Purpose: 

To  compute  the  element  stress  and  stiffness  matrices 
for  a  shear  web  element. 

3.  Input  Arguments: 

X, Y,Z  Node  point  coordinates 

E  Modulus  of  Elasticity 

AMU  Poisson's  Ratio 

RHO  Element  Density 

N5jN6,N7,N8  Element  nodes 
NEL  Element  number 

N1  Degrees  of  freedom/node 

N2  Total  number  of  nodes 

4.  Output  Arguments : 

MM  Unreduced  direction  numbers  for  degrees 

of  freedom  element 

MM(25)  Order  of  stiffness  matrix 

MM(26)=1  Number  of  stress  components  (on  tape  NSS3) 
STfl)  Start  of  stiffness  matrix 

ST(301)  Start  of  stress  matrix 

ARC-1  (On  tape  NSS3)  weight  of  element  with  unit 

design  variable 

5.  Error  Returns:  None 

6.  Calling  Sequence : 

CALL  ELEM2  (X,Y,Z,E,AMU,RH0,N7,N5,N6,N8,N9,N10,N11,N12 
R,MM,NEL,N1,N2) 

7.  Input  Tapes:  None 

8.  Output  Tapes:  NSS3  (Unit  3) 

9.  Scratch  Tapes:  None 

10.  Storage  Required: 

(7246  Bytes)  1811  words 

11.  Subroutine  User:  AONE 

12.  Subroutine  Required:  None 

13.  Remarks: 

A.  This  element  can  be  defined  by  2  or  4  nodes. 

B.  The  stress  matrix  provides  one  stress  component. 


Subroutine  Name:  ELEM3 


Purpose: 

To  compute  the  element  stress  and  stiffness  orthotropic  matrices 
for  a  triangular  plate  element,  with  properties  and  material  angle 
variation. 


Input  Arguments: 


X,  Y.Z 
E 

AMU 

RHO 

N5,N6,  N7 

NEL 

N1 

N2 

EM 

IMAT 

ANGLE 


Node  point  coordinates 
Modulus  of  elasticity 
Poisson's  ratio 
Density 
Element  nodes 
Element  number 
Degrees  of  freedom/node 
Total  number  of  nodes 
Material  properties 
Orthotropic  code 
Material  angle 


Output  Arguments : 


M 

M(25) 
M(26)=3 
ST  (1 ) 
ST(301) 
AREA 


Reduced  direction  number  for  degrees  of  freedom 
of  element 

Order  of  stiffness  matrix 

Number  of  stress  components  (on  tape  NSS3) 

Start  of  stiffness  matrix 
Start  of  stress  matrix 

(On  tape  NSS3)  weight  of  element  with  unit  design 
variable 


Error  Returns:  None 

Calling  Sequence: 

CALL  ELEM3  (X,  Y,  X,  E,  AMU,  RKO,  N7,  N5,  N6,  N8,  N9,  N10,  N1 1,  N12, 
R,  M,  NEL,  Nl,  N2,  GM,  IMAT,  ANGLE) 

Input  Tapes:  None 

Output  Tapes:  NSS3  (Unit  3) 


Scratch  Tapes: 


None 


Storage  Required: 


(4876  Bytes)  1219  words 
AONE 


Subroutine  User: 

Subroutine  Required:  ORTH03 

Remarks: 

Three  components  of  stress  at  the  midpoint  are  provided. 
These  stresses  are  Ox,  (fy,  'fay. 


if; 

!i 

IS 

A 


•\  '  ! 

X  I 

w  i  if  '» 


f  *■ 

it;  1. 

tet% 


1 


1. 

2. 


Subroutine  Name : 
Purpose : 


ELEM4 


This  routine  computes  the  element  stress  and  stiffness 
matrices  for  a  quadrilateral  plate  element. 

Input  Arguments : 


X,Y3Z 

!mu 

RHO 

N5,  N6, 
N7,  N8 
NEL 
N1 
N2 


Node  point  coordinates 
Modulus  of  elasticity 
Poisson's  ratio 
Density 

Element  nodes 

Element  number 
Degrees  of  Freed om/node 
Total  number  of  nodes 


Output  Arguments : 
M 


M  25) 

M(26)=7 

ST(l) 

ST (301^ 
ARGl 

x2/x3 


Unreduced  direction  numbers  for  degrees 
of  freedom  of  element 
Order  of  stiffness  matrix 
Number  of  stress  components  (on  tape  NSS3) 
Start  of  stiffness  matrix 
Start  of  stress  matrix 

(On  tape  NSS3)  weight  of  element  with  unit 
design  variable 

Ratio  (on  tape  NSS3 )  used  in  SCALE  routine 


c; 

>  • 

Error  Returns : 

None 

6 . 

Calling  Sequence : 

i 

CALL  ELEM4  (X,  Y,Z. 
N123R3M,NEL3N13N2) 

s,amu3rho,n7jN5,n63n8,n93nio3nii, 

7.- 

Input  Tapes  : 

None 

,  t 

8. 

Output  Tapes : 

NSS3  (Unit  3) 

9- 

Scratch  Tapes : 

None 

■  % 
i  4 

10. 

Storage  Required: 

(6544  Bytes)  1636  words 

if 

t 

-»! 

11. 

Subroutine  User: 

AONE 

12. 

Subroutine  Required:  MJ]W 

BCB 

'j 

u 

13. 

st 

Remarks:  Seven  components  of  stress  are  defined.  These 

stresses  are: 

^  xi . 


1,  v~iy,  CJxm3  ym3  ^xym  (M-indicates  midpoint) 
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{ 

I 

l 

I 

i 

a 


1. 


2. 


3. 


4. 


5. 


6 . 


7. 


8. 

9. 

10. 


11. 


12. 

13. 


Subroutine  Name :  ELEM5 


Purpose : 

This  subroutine  computes  the  element  stress  and  stiff¬ 
ness  matrices  for  a  beam  (tube)  element. 

Input  Arguments : 


X,Y,Z 

E 

AMU 

RHO 

N5,N6,N7 

NEL 

N1 

N2 


Node  point  coordinates 
Modulus  of  elasticity 
Poisson's  ratio 
Density 
Element  nodes 
Element  number 
Degrees  of  freedom/node 
Total  number  of  nodes 


Output  Arguments: 

M  Unreduced  direction  numbers  for  degrees 

of  freedom  of  element 
M(25)  Order  of  stiffness  matrix 

m(26)=6  No.  of  stress  components  (on  tape  NSS3) 
ST(1)  Start  of  stiffness  matrix 

ST (301)  Start  of  stress  matrix 

ARG1  (On  tape  NSS3)  weight  of  element  with’ unit 

design  variable 


Error  Returns :  None 


Calling  Sequence: 

CALL  ELEM5  (X,Y,Z.E,AMU,RH0,»7,N5,N6,N8,N9,N10,N11, 
N12,R,M,NEL,N1,N2) 

Input  Tapes :  None 

Output  Tapes:  NSS3  (Unit  3) 

Scratch  Tapes :  None 

Storage  Required:  (4466  Bytes)  1116  words 

Subroutine  User:  AONE  ^ 

Subroutine  Required:  3C3 

Remarks : 


Six  stresses  are  available  from  the  stress  matrix. 

These  stresses  are 

crk,  2yz,  3"my1}  crmy2,  Omz1, 


123 


I 
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1.  Subroutine  Name:  ELEM6 

2 .  Purpose : 

This  subroutine  computes  the  element  stress  and  stiff¬ 
ness  matrices  for  an  axial  element  with  a  mid-point 
node , 


3. 


k. 


5. 


Input  Arguments : 


X,Y,Z 

E 

AMU 

RHO 

N5,N6,N7 

NEL 

N1 

N2 


Node  point  coordinates 
Modulus  of  elasticity 
Poisson’s  ratio 
Density 

Element  nodes  N7  (NEL)  is  the  midpoint  node 

Element  number 

Degrees  of  freedom/node 

Total  number  of  nodes 


Output  Arguments : 


M  Unreduced  direction  numbers  for  degrees  of 

freedom  of  element 

Mf25)  Order  of  stiffness  matrix 

M(26)=2  No.  of  stress  components  (on  tape  NSS3) 
ST(l)  Start  of  stiffness  matrix 

ST(301)  Start  of  stress  matrix 

WT  (On  tape  NSS3)  weight  of  element  with 

unit  design  variable 

Error  Returns :  None 


6.  Calling  Sequence: 

CALL  ELEM6  (X, Y, 2, E,AMU, RHO, N7,N5,N6 ,N8,N9,H10, Nil, 
N12,R,M,NEL,  N1,N2) 


-  7 

■1 

7. 

Input  Tapes : 

None 

:  *# 

8. 

Output  Tapes : 

NSS3  (Unit  3) 

*4 

9- 

Scratch  Tapes : 

None 

4 

10. 

Storage  Required: 

(267k  Bytes)  669  words 

!} 

11. 

Subroutine  User: 

AONE 

i; 

1 

12. 

Subroutine  Required:  None 

+  , 


'  if 

i*"  U  t 


•if 


13.  Remartcs : 

A.  Two  stress  components  are  provided  by  the  stress 

matrix.  These  stresses  are:U"x^  andOXg. 

B.  The  coordinates  for  the  mid -point  node  arc  computed 
in  the  routine. 

124 


1.  Subroutine  Name:  ELEM7 

2.  Purpose: 

This  subroutine  computes  she  element  stress  and  stiff¬ 
ness  matrices  for  a  triangular  plate  element  with 
midpoint  nodes . 

3 .  Input  Arguments  : 

X,  Y,Z  Node  point  coordinates 

E  Modulus  of  elasticity 

AMU  Poisson's  ratio 

RHO  Density 

nIJnI’nio!  Element  nodes 

NEL  Element  no . 

N1  Degrees  c*  freedom/node 

N2  Total  no.  of  nodes 

4 .  Output  Arguments : 

M  Unreduced  direction  numbers  for  degrees 

of  freedom  of  element 
Mf25)  Order  of  stiffness  matrix 

M{26)=12  No.  of  stress  components  (on  Tape  NSS3) 
ST(1)  Start  of  stiffness  matrix 

ST (301)  Start  of  stress  matrix 

WT  (On  tape  NSS3)  weight  of  element  with  unit 

design  variable 

5 .  Error  Returns :  None 

6.  Calling  Sequence: 

CALL  ELEM7  (X, Y^.EjAMUjRHOjNTjNSj^NSjNSjNIOjNH, 
N12,R,M,NEL,N1,N2) 

7..  Input  Tapes:  None 

8.  Output  Tapes:  NSS3  (Unit  3) 

9.  Scratch  Tapes :  None 

10.  Storage  Required:  (8022  Bytes)  2006  words 

11.  Subroutine  User:  AONE 

12.  Subroutine  Required:  MINV 

BCB 
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13. 


Remarks : 

A.  Twelve  components  of  stress  are  provided  by  the 
stress  matrix.  These  stresses  are: 

^y-^j  O'Xg.j  >'y2,‘z'xy2,  ixy^3 

<r*B  °5'm> 

where  m  -  denotes  the  centroid  of  the  element. 


B.  Coordinates  of  the  midpoint  nodes  are  computed 
in  the  routine . 


1. 


Subroutine  Name :  ELEM8 

Purpose : 


2. 


3. 


This  subroutine  computes  the  element  stress  and 
stiffness  matrices  for  a  quadrilateral  plate  element 
with  midpoint  nodes . 

Input  Arguments : 


XC,YC,ZC  Node  point  coordinates 
E  "  Modulus  of  elasticity 
AMU  Poisson’s  ratio 


RHO  Density 

N5,N6,N7  ) 

n8,N9,N10V  Element  nodes 
N11,N12  ) 

NEL  Element  no . 

N1  Degrees  of  freedom/node 

N2  Total  no.  of  nodes 


4 .  Output  Arguments : 

M  Unreduced  direction  numbers  for  degrees  of 

M(25)  Order  of  stiffness  matrix 

M(26)=15  No.  of  stress  components  (on  tape  NSS3) 
ST(1)  Start  of  stiffness  matrix 

ST (301)  Start  of  stress  matrix 

WT  (On  tape  NSS3)  weight  of  element  with 

unit  design  variable 

5 .  Error  Returns : 

If  the  ranking  routine  MFGR  outputs  a  rank  not  equal 
to  16,  there  will  be  a  print  out  (3  matrix  error). 
This  will  usually  be  due  to  a  program  input  error. 

6.  Calling  Sequence: 

CALL  ELEM8  (XC3YC3ZC,EJAi.IU3RHO,N7,N5JN6,N8,N9,NlO, 
N11,N12,R,M,NEL,N1,N2) 

7.  Input  Tapes:  None 

8.  Output  Tapes:  NSS3  (Unit  3) 

9.  Scratch  Tapes:  None 

10.  Storage  Required:  (31162  Bytes)  7791  words 

11.  Subroutine  User:  ACNE 
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Subroutine  Required: 


MFCS 

MJNV 

GMERD 

GTPRD 

HOB 

Ksaar -cs : 

A.  Fifteen  stress  components  are  provided  by  the 
stress  matrix-  These  are: 

l'xjk,<r xa,  tx ym. 

where  m  -  denotes  the  element  centroid. 


3.  Coordinates  of  the  midpoints  are  computed  by  the 
subroutine . 
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Subroutine  Name: 


SSCGNT 


Purpose : 

Solves  matrix  equation  A  *  A  transpose  *  x=F  for  X. 
Equations  and  Procedures : 

The  matrix  equation  A  *  A  transpose  X  =  F  is  solved 
for  X  by  solving  the  matrix  equations  A*XX=F 
A  transpose  *  X  =  XX 

where  A  is  a  banded  lower  triangular  matrix  and  X 
and  V  are  column  vectors . 

Procedure: 

fl)  Rows  ISTRT  to  ICALC  of  A  are  read  from  tape  I?. 
(2|  A  call  to  XCALC  routine  computes  XX. 

1 3)  A  call  to  XCALX  routine  computes  X. 

(4)  Steps  1,  2  and  3  are  repeated  for  each  pass. 

Input  Arguments : 

N  Order  of  system 

NPASS  Number  of  computation  passes  necessary 

NEOW  Array  for  control  of  computation  passes 

±z.R  Banding  information  array 

NZEL  Banding  information  array 

A  Storage  array  for  input  matrix 

F  Storage  array  for  input  column 

NTAPE  Input  tape  logical  number 

XX  Wording  storage  array 

NL  No.  of  load  conditions 

Output  Arguments:  X=Output  answer  column  array 

Error  Returns :  None 

Calling  Sequence: 

CALL  SSCGNT  (N,N2ASS,NR0*, L2R,NZ3L,A,F,X,XX,NTAPS,NL) 
Input  Tapes : 

Tape  IT  contains  input  triangular  matrix  A.  A  is 
In  banded  form.  Sac h  row  is  a  separata  record. 

Output  Tapes :  None 

Scratch  Tapes :  None 

:  (2000  Bytes'  500  words 
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Storage  Required 


mmmm 


m 


l 

1. 

Subroutine  Name 

FOMO 

\  j 
[  ; 

2. 

Purpose 

Process 

force  cards 

y 

3. 

Equations  and 
Procedures 

The  cards  are  read  from  file  L5  and 
counted  by  a  counter  KFORCE .  The  input: 

> 

1 

valuer 

are  (Force  l(i),  i  =  1,  2), 

* 

*  * 

■ 

F ,  and 

(F2(i),  i  =  1,  2,  3);  (Force  2  (i )  = 

■  3 

i 

F  -  FnU),  i  -  1,  2,  3). 

1 

i 

! 

4. 

Input  Arguments 

OUT 

=  storage  for  reading  data 

FORCE  1 

=  set  id  number,  grid  point 

number 

i 

] 

FN  (i)  , 

i  =  1 ,  2 ,  3  =  components  of  force 

\ 

L5 

=  input  file  number 

11  ; 

{ 

L6 

=  output  file  number 

b 

[» 

F 

=  scale  factor 

ei 

>; 

si 

si 

5, 

Output  Arguments 

Force  2(i)  for  each  component 

P 

(’ 

i  6. 

Error  Returns 

None 

f 

)* 

7. 

Calling  Sequence 

Call  FOMO  (FORCE  1,  FORCE  2,  L5 ,  NFO , 

1* 

c. 

K  FORCE 

,  OUT,  L6) 

-/ 

sj 

i 

I! 

t\ 

:  ’)•  8’ 

1  ii; 

Subroutine  User 

ZZ 

1.  Subroutine  Name:  GMPRD 

2.  Purpose: 

Multiply  two  general  matrices  to  form  a  resultant 
general  matrix 

3.  Equations  and  Procedures: 

The  M  by  L  matrix  '£»'  is  premultiplied  by  the  N 
by  M  matrix  ’A'  and  the  result  is  stored  in  tne  N  by 
L  matrix  1 R1 . 

4.  Input  Arguments: 

A  Name  of  first  input  matrix 

3  Name  of  second  input  matrix 

N  Number  of  rows  in  A 

M  Number  of  columns  in  A  and  rows  in  3 

L  Number  of  columns  in  B 


5. 

Output  Arguments: 

R  Output  matrix 

6. 

Error  Re  turns : 

None 

7. 

Calling  Sequence-. 

CALL  GMPRD  (A,B,R, 

N,M,L) 

8. 

Input  tapes : 

None 

9- 

Output  Tapes : 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(656  Bytes)  164  words 

12. 

Subroutine  User: 

ELEM8 

13. 

Subroutine  Required:  None 

14. 

Re  mar  Its : 

1.  All  matrices  must  be  stored  as  general  matrices. 

2.  Matrix  R  cannot  be  in  same  location  as  matrix  A. 

3.  Matrix  R  cannot  be  in  same  location  as  matrix  B. 

4.  Number  of  columns  of  matrix  A  must  be  equal  to 
number  of  rows  in  matrix  B. 
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1. 


Subroutine  Name:  GTFRD 


2.  Purpose: 

Premultiply  a  general  matrix  by  the  transpose  of 
another  general  matrix. 

3.  Equations  and  Procedures: 

Matrix  transpose  of  A  is  not  actually  calculated. 
Instead  elements  of  matrix  A  are  taken  columnwise 
ratner  than  row-wise  for  post-multiplication  by 
matrix  B. 

4.  Input  Arguments: 

A  Name  of  first  input  matrix 

B  Name  of  second  input  matrix 

N  Number  of  rows  in  A  and  B 

M  Number  of  columns  in  A  and  rows  in  R 

L  Number  of  columns  in  B  and  R 


5. 

Output  Arguments : 

R  Name  of  output  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  GTPRD  (A,B,R, 

N,M,L) 

8. 

Input  Tapes : 

None 

9- 

Output  Tapes : 

None 

10. 

Scratch  Tapes : 

None 

11. 

Storage  Required: 

(656  Bytes )  164  words 

12. 

Subroutine  User: 

ELEM8 

13. 

Subroutine  Required:  None 

14. 

Remarks : 

1.  Matrix  R  cannot  be  in  the  same  location  as  matrix  A 

2.  Matrix  R  cannot  be  in  the  same  location  as  matrix  B 

3.  All  matrices  must  be  stored  as  general  matrices. 
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Subroutine  Name: 


I  GEN  2 


2.  Purpose: 

Compute  eigenvalues  used  in  buckling  analysis. 

3.  Equations  and  Procedures: 

Sberieir.s  metnod  is  used.  This  consists  of  a  sequence 
of  similarity  transporations  wnich  are  intended  to 
reduce  the  Euclidean  norm  of  the  matrix  under  study 
(A),  so  tnat  it  can  prac ticaliy  speaking,  be  made 
arbitrarily  normal. 

4.  Input  Arguments: 

A  Matrix  wnos.e  eigenvalues  are  desired 

N  Order  of  A-matrix 

YR  Used  to  control  convergence  (min.  IE-7 

smallest  element  of  A) 

ISW=0  Dc  not  form  eigenvectors  (control) 

H)= 6  Dimension  of  'A'  in  calling  routine 

5-  Output  Arguments: 

A  Overlayed  with  matrix  of  eigenvalues  along 

diagonal 

IT  Iteration  counter  (maximum  of  50) 

6.  Error  Returns: 


IIT  If  off  diagonal  elements  of  'diagonalized’ 


7. 

Calling  S  e  que nc  e : 

CALL  IGEN2  (N, YR, 

ISW,  IT,  IIT,  ID,A,AA,T,TI,RES  ) 

8. 

Input  Tapes: 

None 

Q 
-/  • 

Output  Tapes : 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(5700  Bytes)  1425  words 

12. 

Subroutine  User: 

3UCKTB 

13. 

Subroutine  Required:  None 

14. 

Remarks : 

None 
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Subroutine  Name 


INSPC 


Purpose  Stores  data  obtained  from  input  card 

into  working  storage 

Equations  and  SPCINF  =  Out  storage,  SPCTYP 

Procedures 

Input  Argument  SPCTYP  =  Special  code  for  type  of  card 

OUT  =  Storage  for  card  data 

Output  Argument  SPCINF 

NSCARD  =  Counter  for  no.  special  input 
cards 

Error  Returns 

Calling  Sequence  Call  INSPC  (NSCARD,  OUT,  NSPCS,  SPCTYP, 

SPCINF) 

Subroutine  User 


Subroutine  Name 


2.  Purpose 


3.  Equations  and 
Procedures 


4 .  Input  Argument 


INSRT1 

Stores  element  properties  and  material 
numbers  into  working  storage 


F  (I) 

G(I) 

MATNO(I)  = 


If  NA  =  NAST (I)  and 
If  NP  =  NPROP (I) 


NAST 

NPROP 


=  Element  property 
=  Element  property 
=  Material  property  SID 
=  Number  elements 
=  Element  id  number 
=  Property  ~iD 
=  Work  value  of  NAST 
=  Work  value  of  NPROP 


5 .  Output  Argument 


6 .  Error  Return 


None 


7.  Calling  Sequence 


8.  Subroutine  User 


Call  INSRTI  (NAST ,  NA,  NPROP,  NP  MATNO 
NM,  F,  FA,  G,  GA,  NUM) 


1. 


Subroutine  Name: 


KCALC 


2.  Purpose: 

Solve  for  XX  where  A*XK=F  and  A  is  a  banded  lower 
triangular  matrix. 

3-  Equations  and  Procedures: 

De carmine  elements  ISTRT  to  ICALC  of  each  load  vector  XK 
in  the  matrix  equation  A*XX=F. 

F  I  -  1 

XK(1)  =  ~~  XK(I)  =  F(I)  -  t  Att  XK(L)  1 
A11  L  =  1 

This  constitutes  the  first  step  in  equation  solution 

by  Cholesky  or  "square  root"  method. 

Repeat  above  procedure  for  all  load  conditions . 

Input  Arguments: 

ISTRT  Beginning  row  number  of  computation  pass 

ICALC  End  row  number  of  computation  pass 

IZR  Number  zero  elements  in  row  of  reduced  matrix 

NZSL  Cumulative  total  of  nonzero  elements  in  rows 

1  through  i  -  1  of  reduced  matrix 
A  Storage  array  for  input  matrix 

F  Column  vector  array 

N  Order  System 

NL  Number  of  load  conditions 


5. 

Output  Arguments : 

XK  =  output 

vector  array 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

CALL  KCALC  ( ISTRT,  ICALC,  IZR, 

NZEL3A,F,XK,  i 

8, 

Input  Tapes : 

None 

9. 

Output  Tapes : 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(310  3ytes 

)  203  words 

12. 

Subroutine  User: 

ESCONT 

13. 

Subroutines  Requi 

red :  None 
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Subroutine  Name 
Purpose 

Equations  and 
Procedures 

Input  Argument 

Output  Arguments 

Error  Returns 
Calling  Sequence 


KMAT 

Generate  stiffness  matrix  for  membrane 
triangle 

SYMM  =  BINV*  C  BINV * H* AREA 

S  =  Intermediate  work  storage 

BIUV  =  BINV  matrix  (transformation) 

C  =  Elasor  C  matrix 

AREA  =  Area  of  triangle 

H  =  Thickness  of  plate 

KSYMM  =  output  symmetric  stiffness 

matrix 

None 

CALL  KMAT  (S,  BINV,  C,  KSYMM,  AREA,  H) 


Subroutine  User 


0RTH03 


Subroutine  Same :  MFGR 


Purpose : 

Determine  rank  and  linearly  independent  rows  and 
columns  of  a  given  matrix 

Equations  and  Procedures : 

The  standard  Gaussian  elimination  technique  with  complete 
pivoting  is  used.  This  implies  that  the  rows  and 
columns  of  the  given  M  by  N  matrix  A  are  interchanged 
at  each  elimination  step  if  necessary.  The  interchange 
information  is  recorded  in  two  integer  permutation 
vectors  IROVT  and  ICOL. 

The  J  Rcw  *>  of  the  interchanged  matrix  corresponds 
(  Column  ) 

to  the  (^)  ^ow  "l  in  the  original  matrix  A 

^ICOL  (I)°k  Column/ 

Initially  IROTW(J)  =  J  and  ICOL(J)  =  J 

Input  Arguments : 

A  Given  matrix  with  M  rows  and  N  columns 

M  Number  of  rows  of  matrix  A 

N  Number  of  cols  of  matrix  A 

EPS  Test  value  for  zero  affected  by  roundoff  noise 

Output  Arguments : 

IRANK  Resultant  rank  of  given  matrix 

IRCW  Integer  vector  of  dimension  M  containing 

the  subscripts  of  basic  rows  in 
EtOW(l).  .  . IROW(RANK) 

ICOL  Integer  vector  of  dimension  N  containing 

the  subscripts  of  basic  columns  in 
IC0L(1)  UP  TO  ICOL (RANK) 

Error  Returns :  None 

Calling  Sequence: 

CALL  MFGR  (A , M,N, EPS, IRANK, UtOW, ICOL) 

Input  Tapes :  None 

Output  Tapes :  None 

Scratch  Tapes :  None 


1.  Subroutine  Name:  MINV 

2.  Purpose:  Invert  a  matrix 

3.  Equations  and  Procedures: 

Uses  the  standard  Gauss -Jordan  method  in  which  the 
inverted  matrix  is  stored  back  on  itself. 

4.  Input  Arguments : 

A  Matrix  to  be  inverted 

N  Oraer  of  matrix 

D  Determinant  value 

L  Work  vector  of  length  N 

M  Work  vector  of  length  N 

5.  Output  Arguments: 

A  Contains  the  inverted  matrix 

6.  Error  Returns : 

If  D=0,  matrix  is  singular 

7.  Calling  Sequence: 

CALL  MINV  (A,N,D,L,M) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

1C.  Scratch  Tapes'  None 

11.  Storage  Required:  (2C68  Bytes)  517  vrords 

12.  Subroutine  User: 

EL3M4,  ELEM7,  ELEM8 

13.  Subroutine  Required:  None 

14 .  Remarks :  None 
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Subroutine  Name:  MODIFY 


Purpose : 


Evaluate  the  recursion  relationships  for  the  iterative 
determination  of  an  elements  design  variable  when 
a  displacement  constraint  is  violated. 


Equations  and  Procedures : 


^  V~,o 

general 

/_] 

recursion  relatic 

|aq| 

'j  l  F 

U  j  l  JVbJ  |VJ  1 

A 

.1  .  * 

c,_ . 

■  ^1 _ 

u 


k 


Ckjijkj4j 


is 


which  allows  for  all  members  of  a  linked  group  to 
participate  in  the  redesign.  Members  are  considered 
passive  in  the  above  equation  if  the  size  of  that 
member  is  defined  by  other  than  a  displacement  constraint. 
The  largest  area  generated  in  the  above  equation  is  taken 
from  the  combined  set  resulting  from  all  loads  and 
displacement  constraints.  These  areas  are  then  compared 
with  those  based  upon  stresses  or  minimum  sizes  and  the 
larger  values  selected  for  each  member. 


Input  Arguments: 


DELIA 


LIST 


PA?. 


Lower  and  upper  stress  bounds  for  each  element 
Minimum  design  variable 
Upper  and  lower  limits  for  displacement 
constraints 
Displacement:  matrix 

Residual  strength  ratio  output  by  scale 
routine 

Array  containing  number  of  stress 
components  for  each  element 
Element  information  array  originally 
on  NS S3 


Output  Arguments: 

AREA  New  design  variable  for  each  element 


Error  Returns : 


N  one 


Calling  Sequence: 

CALL  MODIFY  (3IGL, SIGU, AIL,  N7,DISFJ,DISPL,  NBDF,  NS 
DELTA, 3IGS ,ARHA , STRESS  LIST, PAR, NBDF2, NE2,AQEA, N5 
N6, X, AREAD, NSTRES, NSS 1) 

Input  Tapes : 

NSS1  (Unit  1)  Element  stiffness  matrices 

Output  Tapes :  None 

Scratch  Tapes:  None 

Storage  Required:  (7516  Bytes)  1879  words 

Subroutine  User:  ATVO 

Subroutine  Required:  None 


Remarics : 


None 


subroutine  Name :  Niwo 


Purucse : 


To  determine  max  core  a  liable  vs  problem  size  and 
set  up  dynamic  storage  subscripts  for  the  initialization 
phase  of  the  program. 

Sanations  and  Procedures: 


Storage  requirements  for  the  arrays  vised  by  the  AONE 
(utilization  routine)  are  computed  from  the  input 
parameters . 

If  the  storage  required  exceeds  MAXCOR  (the  size  of  the 
V  array)  then  ISP.  is  set  to  one,  a  message  is  written 
describing  the  additional  storage  required  and  the 
routine  returns  to  the  main  routine. 

If  the  storage  available  is  adequate,  locations  in  the 
w  array  are  computed  for  the  required  arrays  and 
routine  AC NS  is  called. 


Calling  Argument 


(*  Input) 


*N2 
4  Tarot' 

*USG 

*NAAALL 

*N3CU 

*N3YM 
*NCNS  TA 
*NASYM 

•*KAA 
*NTAPS 
*N?GT 
*NSS  1 ) 
*NSS2 S 
«NSS3i 
•'•NSSh ; 
♦I3UK1* 
*NNZL 


Control  for  degrees  of  freedom  of  a  node 
set  to  3  internally. 

No.  of  nodes  in  problem. 

Number  of  elements  in  problem 
Number  of  components  in  integer  array 
indicating  linked  groups 
Number  of  generalized  constraint  groups 
(NAA)  +  cumulative  length  of  all  groups 
Total  number  of  constrained  degrees  of 
freedom  for  a  symmetric  lead  condition 
Number  of  symmetric  load  conditions 
Number  of  nor.symmetric  load  conditions 
Number  of  anti -symmetric  load  conditions 
Number  of  individual  constraints 
Number  of  generalized  constraints 
Input  tape  i  or  AONE  routine  (Unit  8) 
Print -out  data  set  (Unit  6) 

'Unit  1 

Cutout  data  sets  for  \ 

- Unit  £ 

*fe*.£*i 

Numoer  of  loads  input 
Total  r. ’ember  of  load  conditions 
( NS  YM-rNA  S  YM-rNONS  YM ) 

Number  of  nodes  on  symmetry  plane 
Maximum  number  of  itaratiens* 


Cutout  data  sets  for 
initialization  ohase 


Print  control  for  printing  input 
Integer  array  indicating  number  of  each 
type  element,  e.g.,  ISLI(5)=N  type  5  elements 


NREACT 

*C1INP 

*C2INP 

NRDF 

NRDF2 

NEAA 

NEA2 

NDLI 

NDL2 

NBOUN  2 

*KNLMAX 

*IRST 

*MAXCOR 

*W 

IER 


Integer  array  six  elements  long  indicating 
total  number  of  constraints  e.g.  NREACT(1)= 
number  of  constrained  U  components 
Convergence  control 
Convergence  control 

Number  of  reduced  degrees  of  freedom  for 
symmetric  load  condition 
Number  of  reduced  degrees  of  freedom  for 
antisymmetric  load  condition 
Number  of  generalized  constraint  groups  (NAA) 

+  cumulative  length  of  all  groups 
Same  as  above 

Sum  of  individual  constraints  minus  bounded 
constraints  for  symmetric  load  condition 
Same  as  above  except  for  antisymmetric  load 
condition 

Total  number  of  constrained  degrees  of  freedom 
for  antisymmetric  load  condition 
Max.  length  of  a  generalized  constraint  group 
Calculation  control  which  determines  use  of 
OPDVIR  section 

Size  of  array  available  for  dynamic 
storage 

Array  used  for  dynamic  storage 
Error  control 


5.  Error  Returns: 

lER^O  rot  enough  core  available  to-  complete  initialization 

6.  Calling  Sequence: 

CALL  NEWS (N1,N2 , ITOT,  NSG, NAAALL,  NBOU,  NSYM, NONSYM, NASYM 
NDL, NAA , NTAPE, NPOT , NSS 1 , NS  S2 , NSS  3 , NSS4 , IBUKL , NNZL , NALD , 
KLN2 ,  ITERN,  NTCDS ,  IREST,  IELI,  NREACT,  CULT ,  C2INP .  NRDF ,  NRDF2 , 
NEAA , NEA2, NDLI, NDL2 , NB0UN2 , KNLMAX, IRST, IBGP, MAX COR, W, IER) 

7 .  Input  Tapes :  None 

8.  Output  Tapes:  None 

9.  Scratch  Tapes:  None 

10.  Storage  Required:  (3436  Bytes)  839  words 

11.  Subroutine  User:  MAIN 


12.  Subroutine  Required:  AONE 

13.  Remarks;  None 
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Sub  routine  Name : 


OPINPT 


Purpose : 

To  complete  NTAPE  (Unit  8);  input  to  the  AONE  routine. 
Equations  and.  Procedures  : 

Process  input  needed  by  the  AONE  routine  in  this  order 

COORD 
ELEM 
0 EXTERN 
LIN 1*5 
BOUND 
OLOADS 
ICON 
GCON 

opdvih 

Place  all  information  on  NTAPE  (Unit  8). 

Print  diagnostic  messages  if  necessary. 

Input  Arguments: 

Unit  5  (card  reader) 

Unit  6  (printer) 

Unit  8  (tape) 

Number  of  nodes 

Number  of  individual  constraints 
Number  of  generalized  constraints 
Number  of  input  cards  for  loads 
Unit  1  (tape) 


Info  from  System  section  input 


Output  Arguments : 

NREACT  Array  indicating  number  of  bounded  DOF 

IER  Error  indicator 

Error  Returns :  IER=0 

Calling  Sequence : 

CALL  OPINPT (NP IT, NPOT, NTAPE, NREF, NGRD , NDOFPN , NALD , ITOT 
NIBCP, N1,N2, NSGIN, NDL, NAA , NDLA , NREACT, NNZL, IELI, LRST, 
NAANUM,  NSG,  NSS1,  KNLMAX,  IER,  GRID,  IBOUND,  NBOUND,  INODE, 
NBDF,DIS?U.DISPL, LINKB,LINKN, ILINK,DVIR, IBUCKL, NAM  ,  IBO 
IG1,IG2,IG3,IBGP) 


NPIT 

NPOT 

NTAPE 

N2 

NDL 

NAA 

NNZL 

NSS1 

NREF  ) 

NGRD  ) 

NDOFPN) 

NALD 

ITOT 

NIBCP  ) 


8. 

Input  Tapes : 

None 

9. 

Output  Tapes : 

NTAPE 

10. 

Scratch  Tapes : 

NSS1  (Unit  1)  used 

to  merge  ELEM  and  OEXTRN  sections 

n. 

Storage  Required: 

(15844  Bytes)  3961  words 

12. 

Subroutine  User: 

MAIN 

13. 

Subroutine  Required 

:  None 

14. 

Remarks : 

None 

Subroutine  Name  0PTIM2 


Purpose : 

To  control  the  program  execution. 

Equations  and  Procedures: 

Identify  all  external  files  and  rewind  them. 

Read  REPORT,  TITLE,  SYSTEM  and  OPTIM  sections  to  define 
variables  needed  for  dynamic  storage. 

If  there  are  any  title  cards  read  them  and  put  them  on 
NTAPE  (Uni t  8) . 

Print  message  if  any  sections  are  out  of  order. 

Call  OPINPT  routine  to  read  rest  of  input  sections. 

If  there  is  any  error  (IER^O)  call  exit. 

Call  NEWS  routine  to  perform  dynamic  storage  allocation 
and  call  routines  to  perform  initialization. 

If  there  is  any  error  (IER^O)  call  exit. 

Call  SIZE  routine  to  perform  dynamic  storage  allocation 
and  call  routines  to  perform  problem  calculation. 

If  there  is  any  error  (IER^O)  call  exit. 

Input  Arguments:  None 

Output  Arguments:  None 

Error  Returns:  (IER^O)  call  exit 

Calling  Sequence:  Call  0PTIM2  (WORK,  NP IT) 

Input  Tapes:  WORK  =  Work  storage 

NPIT  =  Input  file  number 

Output  Tapes:  NTAPE  (Unit  8) 

Scratch  Tapes: 

NSS 1  (Unit  1) 

NSS2  (Unit  2)  These  tapes  are  defined  and  rewound 

NSS 3  (Unit  3)  only. 

NSS 4  (Unit  4) 

Storage  Required: 

(35630  Bytes)  8907  words 

This  size  includes  1 H *  work  array  of  8000. 

This  array  is  used  to  dynamically  locate  arrays. 

Subroutine  User:  Hot  applicable 

Subroutine  Required:  OPINPT 

NEWS 

SIZE 


Remarks:  None 


Subroutine  Name : 


FS3D 


Purpose : 


Create  cubic  equation  coefficients  from  stress 
components  for  beam  (tube)  element 


Equations  and  Procedures: 

r3-  (rB+  ryy+  rz2>  ra+(  rxx  ry/+V««.+‘S* 


_<v2_.^-2__/  *)'?'-(  y  "y  ~y  -  'y  'Y  ^  «v  ‘ 
xz  xy  ryz  u  xx^yy^zz  ^xx^yz  *  yy  ^xz 


-  r  r  +  2  2-  ^  T  ) 

zz  'xy  xy  yz  *xz' 


=  0 


or 

,  2 

ry-kl  +  B^-C  =  0 

Input  Arguments : 

SI  -  S6-  Stress  components  (xx,yy,zzpcy1yz jZx) 

XX  Proportioning  factor 

Output  Arguments : 

A  ") 

B  x  Roots  returned  from  cubic 

c  J 

Error  Returns:  None 

Calling  Sequence: 

CALL  PS3D(S1,S2,S3,S4,S5,S6,XX,A,B,C, ) 


Input  Tapes : 

Output  Tapes : 

Scratch  Tapes-: 
Storage  Required: 
Subroutine  User: 
Subroutine  Required: 
Remarks : 


None 

None 

None 

(696  Bytes)  176  words 
PS3R 
C  U3IC 
None 
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Subroutine  Name:  PS3R 

Purpose:  Normalize  Stresses 

Equations  and  Procedures: 

X  =  l/average  stress 

Multiply  all  stresses  by  average  stress 

Input  Arguments:  31  ”*86 

Output  Arguments : 

A,  3j  C  Principal  stresses  of  Beam  (tube)  element 

Error  Returns :  None 

Calling  Sequence : 

CALL  PS3R  (31,  32,33,34,35,36^,3,0) 

Input  Tapes :  None 

Output  Tapes  :  None 

Scratch  Tapes :  None 

Storage  Required:  (66b  Bytes)  166  words 

Subroutine  User:  SCALE 

Subroutine  Required:  PS3D 


Remarks : 


None 


1.  Subroutine  Name  READI 

2.  Purpose  Reads  and  modifies  input  data 

3.  Equations  and  Reads  the  card  data  from  file  J5,  right 

Procedures  adjust  the  data  fields,  counts  each 

data  type,  prints  the  data  and  finally 
stores  the  modified  data  on  file  J6 


4. 

Input  Argument 

LABEL 

= 

array  of  BCD  label  codes 

ILAB 

= 

array  of  label  integers 

ISPECL 

= 

array  of  special  labels 

NILAB 

- 

total  no.  of  labels 

J6 

- 

output  file  number 

NSPECL 

= 

total  no.  special  labels 

5. 

Output  Arguments 

NCARDS 

= 

array  of  card  counters 

L7CASE 

code  to  indicate  that  file 

is  =  7 

6. 

Error  Returns 

None 

7. 

Calling  Sequence 

CALL  READI  (LABEL,  ILAB,  ISPECL,  NCARDS 
NILAB,  NSPECL,  J6,  L7CASE ) 

8. 

Subroutine  User 

SORT 

9. 

Subroutines  Used 

ADJUST 
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1.  Subroutine  Name:  SCALE 

2.  Purpose: 

To  scale  the  stresses,  displacements  anddesign  variables 
of  the  structure.  For  the  in-core  solution 

3.  Equations  and  Procedures: 

A.  Determine  maximum  ratio  BIG  =  actual  stress/allowable 
for  all  elements  and  all  load  conditions. 

For  the  beam  (tube)  element  use  maximum  principle 
stress  from  both  ends  as  allowable. 

For  all  other  elements  use  Von  Mises  reference 
s  tress 

Vbsj  =  ( <5X2-  6/-  <5X  5y+3  rj;y2)* 


B.  Determine  max.  ratio  (BIG)  using  individual  and 
generalized  displacement  constraints 

C.  Scale  element  design  variables  (use  minimum  if 
larger). 

If  buckling  is  specified  use  this  criteria  for 
quadrilateral  selected. 


D. 


& . 


Scale  all  actual  element  stresses 


=  ^act/3lG 


Determine  residual  ratio  (BIGS)  for  each  element 
after  scaling  stresses.  First  use 

minimum  design  var. 

BIGS(I)  =  - : — - — : - 

actual  design  var. 

Then  compare  against  ratio  obtained  as  in  step  B) . 
Use  maximum  for  each  element. 


d.  Input  Arguments: 


LX  No.  of  stresses 

SIGL  Lower  stress  limit  for  each  element 

SIGU  Upper  stress  limit  for  each  element 

ALL  Minimum  design  variable  for  each  element 

IT 7  Not  used 

DISPV  Upper  displacement  limits 

DISPL  Lower  displacement  limits  reduced  deg.  of  free. 

N3DF  Integer  array  locating  RDOF  of  individual 

constraints  symmetric  load  conditions 
NE  Integer  array  locating  RDOF  of  generalized 

constraints  anti-symmetric  load  conditions 
DELIA  Displacement  matrix 
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BIGS 


Residual  ratio  for  each  element  after 
scaling 

STRESS  Stress  matrix 

LIST  Integer  array  giving  no.  of  stress 

components  for  each  element 
PAR  Array  containing  element  weight  factors 

and  the  parameter  x2/x~  for  element  4 
(quadrilateral)  ^ 

NBDF2  Integer  array  locating  RDOF  of  jjidividual 

constraints  for  anti-symmetric  load 
conditions 

NE2  Same  as  NE  except  for  anti -symmetric  load 

condition 

AREA  Not  used 

Lengths  of  sides  of  quadrilaterals 

Constant  Array  -  used  to  control  buckling 
IBK  J  analysis 

5.  Output  Arguments: 

AQEA  Element  design  parameters 

6.  Error  Returns:  None 

7.  Calling  Sequence: 

CALL  SCALE  (LK,SIC-L,SIGU,ALL,N7,DISPU,DISPL,NBDF,NE, 
DELTA,  BIGS ,AREA,  STRESS ,  LIST,  PAR, NBDF2,NE2, AQEA, AL,  BL, 
IBUKL, I3K) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes :  None 

11.  Storage  Required:  (9110  Bytes)  2277  words 

12.  Subroutine  User:  ATWD 

13.  Subroutine  Required: 


AL  ) 
BL  J 
I3UKL 


PS3R 

BUCKTB 


Subroutine  Name : 


SCALE1 


Purpose : 

To  scale  'the  stresses  >  displacements  and  design 
variables  of  the  structure  for  an  "out-of -core -solution" . 

Procedure : 

The  procedure  followed  here  is  identical  to  the  SCALE 
routine  exceot  the  stress  matrix  is  on  data  set  NSS3 
(one  column/record). 

A.  Determine  max.  ratio  (3IG)=  actual  stresses/allowable 
stress  for  all  elements  and  all  load  conditions. 

Per  the  beam  (tube)  element  use  matrix  principle 
stress  from  both  ends  as  allowable 
?or  all  other  elmenss  use  Von  Misses  Reference 
stress , 

3.  Determine  max.  ratio  (BIG)  using  individual  and 
generalized  displacement  constraints . 

C.  Scale  element  design  variables  (use  specified  minimum 

(if  larger). 

If  buckling  is  specified  use  this  criteria  for 
quadrilaterals  selected. 

D.  Scale  all  actual  element  stresses. 

S  =  Sact/BIO 

E.  Determine  residual  ratio  (BIGS )  for  each  element 
after  scaling  stresses 

Input  Arguments : 

See  SCALE  write-up  for  more  complete  listing  of  Arguments. 
LK  No.  of  stresses 

SIGL  Lower  stress  limit-  for  each  element 

SIGH  Upper  stress  limit  foi  each  element 

ALL  Minimum  design  variable  for  each  element 

DI3FU  Upper  displacements  limits 

DISPL  Lower  displacement  limits 

DELTA  Displacement  matrix 

NSS3  Data  set  containing  stress  matrix 

Output  Arguments: 

AREA  Element  design  parameters 

Error  Returns:  None 
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Calling  Sequence: 


CALL  SCALS1  (LK,SIGL, SIGU,ALL,  N7,EISPU,DISPL, NEDF, NS, 
DELTA ,  BIGS ,  AREA ,  STRESS ,  LIST,  PAR,  NBDF2 ,  NE2 ,  AQEA ,  AL,  BL, 
IBUKL,  IBK,  NSS  3 ,  BIG ) 


Input  Tapes : 


NSS 3  Stress 

matrix 

Output  Tapes: 

None 

Scratch  Tapes : 

None 

Storage  Required: 

(8836 

Subroutine  User: 

ATW01 

(1  column/record) 


Bytes)  2209  words 


Subroutine  Required: 

PS3R 

BUCKTB 


S ubr  cu ti ne  Name : 


SIZE 


Z  sA  —  » 

Tc  de -ermine  max.  cere  available  vs  problem  size  and 
set  up  dynamic  s corage  subscripts  for  toe  calculation 
pnase  cf  the  program  in-ccre  solution. 

Scua cions  and  Procedures: 


Storage  requirements  for  the  arrays  used  by  the  ATWO 
(calculation  routine)  are  computed  from  the  input 


parameters . 

If  the  storage  required  exceeds  MAXCOR  (the  size  of  the 
X  array)  then  a  message  is  written  describing  the 
additional  storage  required  and  the  program  calls  the 
SIZE!  routine  for  an  attempt  at  an  out-of-core 
solution. 


If  the  storage  available  is  adequate,  locations  in  the 
' x 1  array  are  computed  for  one  required  arrays  and 
routine  ATXC  is  called. 


Input  Arguments: 

Total  number  of  elements 
No.  of  reduced  degrees  of  freedom  for 
symmetric  load  conditions 
No.  of  symmetric  lead  conditions 
No.  of  antisymmetric  load  conditions 
No.  of  ncnsymr.etric  load  conditions 
Printer  data  set 
Data  sets 

Print  controls 

Size  of  *Wf  array 
Array  for  dynamic  storage 

Output  A r gume nts : 

IZR=C  indicates  error  (problem  exceeds  storage  reserved) 
Error  Returns: 


N3LSM 

NRDF 

NSYM 


NCNSYM 

NPCT 

NSSl-NSSp 

tt5  •:  —dc  ) 

IPRINT  ) 

MAXCOR 

W-WCRX 


If  MAXCOR  is  exceeded  then  the  ATXO  routine  is  not 
called  and  tr.e  program  writes  storage  required.  Then 
a  call  is  made  to  the  3IZE1  routine. 

Calling  Sequence : 


CALL  SIZE  (NSLEM,NRDF,Nmi,NAA,NRDF2,NDL2,NSYM,NA5YM, 
NCNSYM.  ISLE,  NSC-,  NEAA ,  NSA2,  NDL,  NBOU,  NB0UN2, NPOT, NSSI, 
NSS2, NSS3 ,NSS4 , NSSp,  I?1,I?2, IRST, ITERN, IPRINT, IBUKL, 
MAXCOR . W, IER ) 
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None 


Input  'Tapes : 

Output  Tapes :  None 

Scratch  Tapes:  None 

Storage  Required:  (2780  Bytes)  695  words 

Subroutine  User:  MAIN 

Subroutine  Required: 

A  "VO 
SI2.TI 


Remarics : 


None 


Subroutine  Name:  SIZZi 


2.  Furpose: 

To  determine  max.  core  available  vs  problem  size  and 
set  up  dynamic  storage  subscripts  for  the  calculation 
phase  of  tne  program  cut -cf -core  solution. 

3-  Equations  and  Procedures: 

Storage  requirements  for  the  arrays  used  by  the  ATW01 
(out-of-core)  calculation  routine  are  computed  from 
the  input  parameters . 

If  the  storage  required  exceeds  MAXCOR  (the  size  of  the 
W  array)  then  a  message  is  written  describing  the 
additional  storage  required  and  tne  program  halts  after 
a  return  to  the  MAUI  routine. 

If  the  storage  available  is  adequate,  locations  in 
the  ‘W1  array  are  computed  for  the  required  arrays 
and  routine  AftfCI  is  called. 


4 .  Input  Arguments  : 


NELSM 

NRDF 

NS  YM 
NASYM 
NCNSYM 
NPOT 

NSS1-NSS5 
IP1,I?2,) 
xrRINT  ) 
MAXCOR 
W 


Total  number  of  elements 
No.  of  reduced  degrees  of  freedom  for 
symmetric  load  conditions 
No.  of  symmetric  load  conditions 
No.  of  antisyrsaetric  load  conditions 
No.  of  nonsymmetric  load  conditions 
Printer  data  set 
Data  sets 

Print  controls 

Size  of  'W*.  array 

Work  array  for  dynamic  storage 


5 .  Output  Arguments : 


IER=0  indicates  error  (problem  exceeds  storage  reserved) 
Error  Returns : 


If  MAXCOR  is  exceeded  then  the  A  TWO  routine  is  r-  t 
called  and  the  program  writes  storage  required.  A 
return  to  the  'MAIN1  routine  then  halts  the  program. 


Calling  Sequence: 

CALL  SIZZI  f  NZLEM,  NRDF ,  NDL1 ,  NAA ,  NRDF2 ,  NDL2 ,  NSYM,  NAS  YM, 
NCNSYM,  IZLI,‘NSG.  NEAA ,  NEA2,NDL,  N30U ,  NB0UN2  ,  NPOT ,  NSS 1 . 
NSS2 , NSS  3 , NSSh , NSS 5 , 1? 1 , IP2 . IRST, ITZSN , IPRINT , IBUXL, 
MAXCOR, w, IZR) 


3.  Innut  Tates : 


None 
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9- 

Output  Tapes : 

None 

10. 

Scratch  Tapes: 

None 

li. 

Storage  Required: 

(2728  Bytes)  682 

12. 

Subroutine  User: 

SIZE 

13. 

Subroutine  Required: 

ATW01 

I 


9 


words 


i 
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1 .  Subroutine  Name 


SMAT 


Generate  stress  element  matrix  for 
membrane  triangular  plate  with  ortotropic 
properties 

SMAK  =  T*C*BINV  where  T  is  generated  in 
SMAT,  C  and  BINV  are  input 

C  =6x6  matrix 
BINV  =6x6  matrix 

5.  Output  Arguments  SMAK  =3x6  stress  matrix 

6.  Error  Returns  None 

7.  Calling  Sequence  CALL  SMAT  (C,  SMAK,  BINV) 

8.  Subroutine  User  ORTHOS 


2.  Purpose 


3.  Equations  and 
Procedures 


4 .  Input  Arguments 


Subroutine  NAME 


Sort 


Purpose 

Equations  and 
Procedures 

Input  Arguments 

Output  Arguments 

Error  Returns 
Calling  Sequence 

Subroutine  User 
Subroutines  Used 


Sort  and  count  data  based  on  LABEL 
information. 


Using  READI,  the  data  is  read  and  counted. 
The  final  counters  are  then  modified. 

L4  =  file  number  for  storage  of  sorted 
input  deck 


*  MEL  =  total  no. 
MGR  =  total  no. 

MMT 

MCGCON  = 

«  MIC0N8  = 

MILINKS  = 

MFO  =  "  " 

MMO 

MLOADS  =  ”  " 

None 

Call  Sort  (MEL,  MGR 
MICON8 ,  MILINKS,  MFO 
MSPCS ,  L4 ,  L7CASE) 


elements 
grid  points 
materials 

generalized  constraints 

INdividual  Constraints 

Links 

Forces 

Moments 

Loads 


MMT,  MOGCON, 

MMO,  MLOADS,  MGROUF , 


Main  Program 
READ  I 


1. 

Subroutine  NAME 

SPCSUB 

2. 

Purpose 

Process  SPC  (single  point  constraint) 

Cards 

Boundary  information  is  processed  as 
read  in  OUT  (I)  and  NOUT  (I).  This 
information  is  interpreted  and  stored 
in  LBOUND. 

4.  Input  Arguments  NOUT  =  input  data  storage 

OUT  =  input  data  storage 

NSPC  =  no.  SPC  cards 

KWORP  -  work  storage 

NUM  =  work  storage 

NKIND  =  type  of  boundary  information 

available 

NGR  =  -total  no.  grid  points 

5.  Output  Arguments  LBOUND  =  boundary  array  information 

NBCARD  =  counter  of  boundary  information 

6.  Error  Returns  None 

7.  Calling  Sequence  Call  SPCSUB  (NOUT,  OUT,  NSPC,  KWORD, 

NUM,  LBOUND,  NBCARD,  NKIND,  NGR) 

8.  Subroutine  User  ZZ 

9.  Subroutines  Used  XTRAK 


3.  Equations  and 
Procedures 
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Subroutine  Name: 


TCONTX 


Purpose : 

This  routine  controls  tape  flow  for  the  triangular- 
ization  routine. 

Equations  and  Procedures : 

1.  Controls  for  setting  up  computation  passes  are  computed 
in  ICALC  and  ISTRT . 

2.  A  portion  of  the  input  matrix  A  is  read  in  from 
Tape  I3=MTAPE. 

3.  This  information  is  given  to  the  routine  TTRI  which 
actually  performs  the  triangularization  for  row 
numbers  ISTRT  to  ICALC. 

4 .  This  triangularized  output  portion  in  A  is  stored 
on  tape  NTAPE  =  17. 

5.  Computation  is  repeated  for  each  portion  of  the 
matrix  until  all  rows  are  completed. 

Input  Arguments : 

N  order  of  system  to  be  handled 

IZR  banding  information  array 

NZEL  banding  information  array 

A  storage  array  for  input  row.  of  banded 

matrix  which  is  read  by  routine 
NTOTAL  total  number  of  words  which  can  be 

considered  as  a  "^ull-core" 

ATRI  intermediate  storage  array  equals 

length  of  maximum  order 
MTAPE  input  tape  logical  number 

NTAPE  output  tape  logical  number 

Output  Arguments : 

IERROR  error  indication  value 

WS  accumulative  determinant 

Error  Returns: 

TERROR  not  =  0  if  WS  is  returned  from  TTRI  as  less 
than  zero 

Calling  Sequence: 

CALL  TCONTX (N, IZR, NZEL, A , NTOTAL, ATRI, MTAPE, NTAPE, 

IERROR, WS) 


Input  Tapes : 

MTAPE  *  13  =  input  matrix  A  in  banded  row  form. 
Each  row  equals  1  record. 

Output  Tapes : 

NTAPE  a  17  a  triangular ized  matrix  T  in  banded  row 
form.  Each  row  equals  1  record. 

Scratch  Tapes :  None 

Storage  Required:  (1592  Bytes)  398  words 

Subroutine  User:  ATW01 

Subroutines  Required:  TTRI 


1 .  Subroutine  NAME 

2 .  Purpose 

3.  Equations  and 

Procedures 

4 .  Input  Arguments 

5 .  Output  Arguments 

6 .  Error  Returns 

7.  Calling  Sequence 

8 .  Subroutine  User 


TMftT 

Generate  transformation  matrix  for 
triangular  plate  orthotropic 
material  angle. 

T  =  fcos2^  Sin2-*-  SineCose- 

Sin2-e  Cos2-*-  -Sin^Cos*- 

2 

-2  .Sin*Cos*  2  SineCos*  Cos  Sin 
THETA  =  material  angle 
T  =  3x3  matrix 

None 

CALL  THAT  (THETA,  T) 

ORTHO 3 


9.  Subroutines  Used  None 
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1. 


Subroutine  Name: 


TTRI 


2.  Purpose: 

To  triangularize  rows  ISTRT  to  ICALC  of  a  bonded 
matrix  A . 

3.  Equations  and  Procedures: 

1.  This  routine  triangularizes  rows  ISTRT  to  ICALC 
of  a  banded  matrix  A  where  rows  1  to  ISTRT-1 

of  the  A  matrix  (already  triangularized)  are  on 
tape  NTAPE. 

2.  If  ISTRT  =1,  then  NTAPE  and  work  storage  ATRI 
are  not  used  since  A  is  assumed  to  be  in  core. 

3.  Procedure: 

Using  Cholesky  technique,  the  off  diagonal 
terms  of  the  portion  in  core  are  triangularized . 
Off  diagonals  are  then  computed.  Output  is  stored 
in  array  A . 

4.  Cholesky  equations: 

(1)  slx  =  (3)  ^Li  =(aii  ’4 S 6£i)^  i 


St0  ■  sjf 


s.  ,=a.  ■ 
ij  i-j 


S£iSJlj 


j  >i 


4.  Input  Arguments : 

ISTRT  Beginning  row  of  triangularized  portion 

ICALC  End  row  of  triangularized  portion 

IZR  3andi_ng  Information  array 

NZEL  Banding  information  array 

NTAPE  logical  tape  number  of  input  tape .NTAPE=I7 

A  Storage  array  for  input  A  and  also  output 

array 

ATRI  Working  storage  array 

5.  Output  Arguments:  A=output  array 

6.  Error  Returns: 


IERROR  =  I  =  row  numbers  such  that  WS  =  s  is  not 
greater  than  zero.  a 
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7.  Calling  Sequence: 

CALL  TTRI  (ISTRT, ICALC, IZR,NZEL, NTAPE,A,ATRI, IERROR,WS ) 

8.  Input  Tapes: 

NTAPE  =  tape  which  contains  already  triangularized  rows 
of  matrix. 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (1844  Bytes)  46l  words 

12.  Subroutine  User:  _  TCONTX 

13.  Subroutines  Required:  None 
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Subroutine  NAME 
Purpose 

Equations  and 
Procedures 

Input  Arguments 

Output  Arguments 
Error  Returns 
Calling  Sequence 
Subroutine  User 


WRITEL 

Tests  character  of  element  connection 
card  and  writes  element  information 
into  file. 

If  and  are  0,  C2  and  are 
stored. 

If  C2  &  C4  are  0,  C-j  &  C3  are  stored. 

L7  =  Output  file  number 
IGR  =  Grid  point  number 

M  =  Position  of  A  array  to  be  restored 
C]_,  C2 ,  C3,  C4  =  Input  codes 

IGR  and  A  array  are  stored  on  file  L  7 

None 

CALL  WRITEL  (Cj,  C2 ,  C3,  C4,  M,  IGR,  L) 
ZZ 


Subroutines  Used 


None 


1. 


Subroutine  Name: 


XCALK 


2.  Purpose: 

Solves  for  X  where  A*X  =  XK  and  A  is  upper  triangular 
matrix. 

3.  Equations  and  Procedures: 

Determine  ISTRT  to  ICALC  of  each  load  vector  X  in  the 
matrix  equation  A*X  =  XK  where  A  is  upper  triangular 
matrix  and  X  and  XK  are  column  vectors. 


X(N) 


X(I)  =  XK(I) 

Ann 


N 

£  A_Xr  I<N 
tIt^i  IL  L 


This  constitutes  the  second  part  in  calculating  an 
equation  solution  by  Cholesfcy  or  "square  root'*  method. 
Repeat  above  procedure  for  all  lead  conditions. 

4 .  Input  Arguments : 

N  Order  of  system 

ISTRT  Beginning  row  number  of  c  deputation  pass 

ICALC  End  row  number  of  computation  pass 

IZR  Number  zero  elements  in  row  of  reduced  matrix. 

NZEL  Cumulative  total  of  nonzero  elements  from 

row  1  through  1  -  iof  reduced  matrix 
A  Storage  array  for  matrix 

XK  Column  vector  array 

NL  No.  of  load  conditions 


5. 

Output  Arguments: 

X  Output  vector  array 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

CALL  XCALK  (N, ISTRT, 

ICALC ,  IZR,  NZEL, A ,  XK, X,  NL) 

8. 

Input  Tapes  : 

None 

9- 

Output  Tapes : 

None 

10. 

Scratch  Tapes : 

None 

li. 

Storage  Required: 

(1116  Bytes)  279  words 

12. 

Subroutine  User: 

ESCONT 

13- 

Subroutines  Required 

:  None 
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1 .  Subroutine  NAME 


XTRAK 


j 

i 


s 

I 

« 

i 


! 


i 


2 .  Purpose 


Interpret  degree  of  freedom  information 


3.  Equations  and 
Procedures 


4 .  Input  Arguments 


Interprets  NWORD  and  breaks  this  down 
into  6  individual  components.  These 
components  are  then  stored  in  KWORD 
array. 

NWORD  =  No.  components  input  word 
NP  =  Control  word 


5 .  Output  Arguments 


KWORD  =  Output  data  array 

NUM  =  Total  no.  DOF  recognized 


6.  Error  Returns  None 


7 .  Calling  Sequence 


CALL  XTRAK  (NWORD,  KWORD,  NUM,  NP) 


8.  Subroutine  User  SPCSUB 


9.  Subroutines  Used  None 
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Subroutine  NAME 
Purpose 

Equations  and 
Procedures 

Input  Arguments 

Output  Arguments 
Error  Returns 
Calling  Sequence 
Subroutine  Oser 


XTRAK2 

Interpret  degree  of  freedom 
information  for  generalized  constraint 
data. 

Checks  on  information  supplied  in 
NP  and  stores  codes  into  LGCON  array. 

NP  =  Input  word  to  be  interpreted 
NM  =  2nd  input  word  to  be  interpreted 
NODE  =  Node  point  number 

LGCON  =  Generalized  constraint  array 

None 

CALL  XTRAK2  (LGCON.  NODE,  NP,  NM) 


Subroutines  Used 


XTRAK 


1 .  Subroutine  NAME 


ZZ 


2  -  Purpose 

3.  Equations  and 
Procedures 

4 -  Input  Arguments 

5 .  Output  Arguments 

6 .  Error  Returns 

7 .  Calling  Sequence 


8 .  Subroutine  User 

9.  Subroutines  Used 


Generates  OPTIM  data  which  is 
input  by  NASTRAN  format  input  cards . 

Each  card  is  r^ad  and  interpreted 
based  on  content  and  use  in  the  OPTIM 
program. 

L  5  -  Input  file  tape  number 
L  7  =  Output  file  tape  number 

All  of  the  grid  point,  boundary 
condition,  element,  material  property 
load,  constraint  and  buckling 
information  arrays  needed  by  OPTIM. 


CALL  ZZ  (NAST,  NOPT,  MATNO,  NPRQP, 
NBUQK,  NNODES ,  NOID,  REF,  OPDVIR, 
LBOUND,  COOR,  AMAT,  MID,  EYEC,  LINKS, 
NEL,  NGR,  NMT,  NICON8,  NLINKS,  NOGCON 
NP,  FA,  NAP,  GA,  NFO,  NHG,  NLOADS , 
NSPCS,  IGRID,  SPCINF,  GCOND,  FORCES, 
ANGLE,  FORCEl ,  FORCE 2 ,  MOMNTl,  MOMNT2 
OPLOAD,  GROUP,  L5,  L7). 

Main  program 

Hone 
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SECTION  5 


PROBLEM  SIZE  LIMITATIONS 


The  key  limitations  are  brought  about  by  the  total  number  of  elements, 
number  of  load  conditions  and  the  number  of  degrees  of  freedom  in  the 
problem. 

Due  to  dynamic  storage  allocation  techniques  the  problem  size  can 
be  controlled  by  changing  two  cards  in  the  'MAIN'  routine.  Both,  the 
dimension  of  the  'WORK'  array  (card  MAIN0050)  and  the  size  of  the  variable 
'NWORK'  (card  MAIN  0290)  must  be  equal.  The  delivered  size  is  20,  000 
words,  but  may  be  adjusted  to  your  system. 

If  there  is  insufficient  storage  space  defined  for  a  problem  the 
program  will  print  a  message  indicating  the  amount  of  storage  required 
for  the  problem  and  the  amount  of  storage  reserved  by  the  above  mentioned 
cards  in  the  'MAIN'  routine.  To  execute  the  problem  modify  thf  ’'mension 
and  'NWORK'  variable  to  be  what  the  problem  needs  or  reduce  the  number 
of  elements,  number  of  load  conditions  and/or  number  of  degrees  of  freedom 
indicated  on  the  input  sheet. 
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SECTION  6 


STORAGE  ALLOCATIONS 

This  section  includes  two  diagrams  indicating  the  manner  of 
arranging  the  stress  and  displacement  matrices  to  handle  all  load 
conditions. 

The  key  to  handling  non- symmetric  loads  and  displacements  is  in 
defining  the  nonsymmetric  load  as  a  symmetric  and  an  anti. -symmetric 
load. 

The  definition  of  a  number  of  key  parameters  used  in  the  diagrams 

follows: 

NDL1  -  Total  no.  of  individual  constraints  minus  any 

constraints  that  may  be  bounded  during  execution 
of  the  problem.  For  symmetric  load  conditions. 

NDL2  -  Same  as  above  except  for  anti -symmetric  load 
conditions. 

NAA  -  No.  of  generaliz  id  constraints,  all  are  assumed 
unbounded. 

NRDF  -  No.  of  reduced  degrees  of  freedom  (i.e.,  all  bounded 
degrees  of  freedom  removed)  for  symmetric  load 
condition. 

NRDF2  -  Same  as  above  except  for  anti -symmetric  load 
condi  lion. 


See  Figures  4- and  o. 
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STRESS  MATRIX  STORAGE  ALLOCATION 


DISPLACEMENT  MATRIX 


KX0AD1 


KL0AD2 


LOAD1 


Lx.' AD? 


NBL1 

NSYM  NONSYM  +NAA 


NDL2 

NASYM  NONSYM  +NAA 


r- 


Symmetric 

Load 

Conditions 


>  NRDF2 
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